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Preface 


The  purpose  of  this  study  was  to  explore  the  structural  and  aeroelastic  behavior 
of  a  standard  swept  vertical  tail.  An  optimization  algorithm  was  used  to  minimize  the 
weight  of  a  design.  Several  different  designs  were  considered  by  looking  at  different 
cases  which  varied  the  design  load,  design  variable  de^nition,  and  design  constraints. 

It  is  hoped  that  the  qualitative  results  of  this  thesis  can  be  easily  applied  to  other 
lifting  surface  designs. 

The  software  that  was  used  for  this  design,  MSC/NASTRAN,  is  based  on 
theory  that  was  not  fully  explained.  Therefore,  a  major  thrust  of  work  for  this  thesis 
was  to  research  the  theoretical  basis  of  this  program  and  understand  and  explain  the 
restrictions  of  the  structural  and  aerodynamic  theory.  The  optimization  algorithm  was 
based  on  a  functional  approach  that  needed  explanation  so  that  the  results  from  the 
iterative  scheme  could  be  closely  monitored. 

I  wish  to  give  special  thanks  to  my  thesis  advisor.  Dr  Peter  Torvik,  and  my 
thesis  sponsor.  Dr  Jim  Olsen,  for  their  guidance  in  times  of  need.  A  word  of  thanks  is 
also  owed  to  Col  Ron  Bagley  for  his  assistance  and  pat  on  the  back  when  I  was 
making  no  headway  on  my  thesis. 


David  G  Miller 


ii 


Table  of  Contents 


Page 


Preface 

ii 

List  of  Figures 

V 

List  of  Tables 

X 

List  of  Symbols 

xi 

Abstract 

xiv 

I.  Problem  Statement  and  Background 

1-1 

Problem  Statement 

1-1 

Introduction 

1-2 

Prediction  of  Aerodynamic  Behavior 

1-3 

Prediction  of  Structural  Behavior 

1-5 

Optimal  Theory  Applied  to  Aeroelastic  Design 

1-6 

Closing 

1-7 

n.  Theory 

2-1 

Introduction 

2-1 

Aerodynamic  Theory 

2-2 

Structural  Analysis 

2-10 

Aeroelastic  Analysis 

2-14 

Structural  Optimization 

2-23 

m.  Vertical  Tail  Model 

3-1 

Introduction 

3-1 

Structural  Model 

3-2 

Optimization  Model 

3-8 

Aerodynamic  Model 

3-16 

Comments  on  MSC/NASTRAN 

3-18 

iii 


IV.  Results 


4-1 


Introduction  4-1 

Effect  of  Design  Variable  Linking  on  Objective  Function  4-1 

Effect  of  Design  Load  Condition  on  Minimum  Weight  4-6 

Effect  of  Constraints  on  Final  Weight  4-10 

The  Static  Aeroelastic  Properties  of  Constrained  Designs  4-16 

Effect  of  Design  Constraints  on  Flutter  Speed  4-19 

V.  Discussion  and  Recommendations  5-1 

Introduction  5-1 

Optimization  5-2 

Aeroelastic  Analysis  5-8 

Closing  Comments  5-9 

Bibliography  bib-1 

Appendix  A:  Modal  Shapes  of  Five  Design  Cases  a-1 

Appendix  B:  MSC/NASTRAN  Input  Bulk  Data  Deck  b-1 

Vita  v-1 


iv 


List  of  Figures 


Figure  Page 

2-1.  Orientation  of  aerodynamic  velocity  vectors  used  in  equation 

formulation .  2-2 

2-2.  Definition  of  kernel  function  coordinate  system .  2-6 

2-3.  Illustration  of  doublet-lattice  mesh .  2-7 

2-4.  Two  1-D  examples  of  interpolation  and  extrapolation  curves .  2-17 

2-5.  Example  optimization  design  space .  2-26 

2- 6.  Graphical  representation  of  0F,  3G,  and  (S) .  2-28 

3- 1.  Vertical  tail  dimensions .  3-2 

3-2.  Equivalent  transfer  of  aerodynamic  loads  from  tail  surface  to  tail  box.  3-4 

3-3.  Description  of  the  six  structural  member  groups .  3-6 

3-4.  Element  types  used  in  design  model .  3-7 

3-5.  Illustration  of  six  and  twenty  one  design  variable  linking  scheme  :  Spar 

(bar)  element .  3-10 

3-6.  Illustration  of  six  and  twenty  one  design  variable  linking  scheme  :  Spar 

panel  element .  3-10 

3-7.  Illustration  of  six  and  twenty  one  design  variable  linking  scheme  :  Rib 

(bar)  element .  3-11 

3-8.  Elustration  of  six  and  twenty  one  design  variable  linking  scheme  :  Rib 

panel  element .  3-11 

3-9.  Illustration  of  six  and  twenty  one  design  variable  linking  scheme  :  Post 

(bar)  element .  3-12 

3-10.  Illustration  of  six  and  twenty  one  design  variable  linking  scheme  :  Skin 

panel  element .  3-12 

3-11.  Aerodynamic  panel  mesh .  3-17 


V 


3- 12.  The  three  NASTRAN  modules  utilized  for  the  vertical  tail  analysis.  .  .  3-19 

4- 1.  Optimal  design  variable  size  -  cross-sectional  area  of  spar  (bar) 

elements .  4-3 

4-2.  Optimal  design  variable  size  -  thickness  of  spar  panel  elements.  .  .  .  4-3 

4-3.  Optimal  design  variable  size  -  cross-sectional  area  of  rib  (bar) 

elements .  4-4 

4-4.  Optimal  design  variable  size  -  thickness  of  rib  panel  elements  ....  4-4 

4-5.  Optimal  design  variable  size  -  cross-sectional  area  of  post  (bar) 

elements .  4-5 

4-6.  Optimal  design  variable  size  -  thickness  of  skin  panel  elements.  .  .  .  4-5 

4-7.  Von  Mises  stress  values  -  Mid  spar  panel  elements .  4-6 

4-8.  Von  Mises  stress  values  -  forward  rib  panel  elements .  4-6 

4-9.  Von  Mises  stress  values  -  forward  skin  panel  elements .  4-7 

4-10.  Optimal  structural  weight  versus  design  angle  of  attack .  4-9 

4-11.  Optimal  weight  of  design  compared  to  applied  aerodynamic  load .  4-9 

4-12.  Design  variable  sizes:  cross-sectional  area  of  spars .  4-11 

4-13.  Design  variable  sizes:  thickness  of  spar  panels .  4-11 

4-14.  Design  variable  sizes:  cross-sectional  area  of  ribs .  4-12 

4-15.  Design  variable  sizes:  thickness  of  rib  panels .  4-12 

4-16.  Design  variable  sizes:  cross-sectional  area  of  posts .  4-13 

4-17.  Design  variable  sizes:  thickness  of  skin  panels .  4-13 

4-18.  Comparison  of  stress  values  in  identical  skin  panel  elements .  4-14 

4-19.  Displacements  at  nodal  points  along  the  three  spars  due  to  rigid  tail 

vi 


loads,  (stress  constrained  case) .  4-18 

4-20.  Displacements  at  nodal  points  along  the  three  spars  due  to 

aeroelastically  corrected  tail  loads,  (stress  constrained  case).  ...  4-18 

4-21.  Displacements  at  nodal  points  along  the  three  spars  due  to  rigid  tail 

loads,  (minimum  gauge  case) .  4-19 

4-22.  Displacements  at  nodal  points  along  the  three  spars  due  to 

aeroelastically  corrected  tail  loads,  (minimum  gauge  case) .  4-19 

4-23.  Displacements  at  nodal  points  along  the  three  spars  due  to  rigid  tail 

loads,  (displacement  constraint  case) .  4-20 

4-24.  Displacements  at  nodal  points  along  the  three  spars  due  to 

aeroelastically  corrected  tail  loads,  (displacement  constraint  case)  .  4-20 

4-25.  Comparison  of  Von  Mises  stress  values  in  mid  spar  panels  with  rigid 

aerodynamic  loads  and  aeroelastically  corrected  loads  applied .  4-21 

4-26.  Comparison  of  Von  Mises  stress  values  in  forward  rib  panels  with  rigid 

aerodynamic  loads  and  aeroelastically  corrected  loads  applied .  4-21 

4-27.  Normal  distribution  of  concentrated  mass .  4-23 

4-28.  Mass  forward  distribution  of  concentrated  mass .  4-23 

4-29.  Mass  aft  distribution  of  concentrated  mass .  4-24 

4-30.  Velocity  -  artificial  damping  (V-g)  curves  for  the  first  4  modes  in  the 

stress  constraint  case .  4-27 

4-31.  Reduced  frequency  -  frequency  (k-OJ)  curves  for  the  first  4  modes  in  the 

stress  constrained  case .  4-27 

4-32.  Velocity  -  artificial  damping  (V-g)  curves  for  the  first  4  modes  in  the 

minimum  gauge  constraint  case .  2-28 

4-33.  Reduced  frequency  -  frequency  (k-GI)  curves  for  the  first  4  modes  in  the 

minimum  gauge  constraint  case .  4-28 

4-34.  Velocity  -  artificial  damping  (V-g)  curves  for  the  first  4  modes  in  the 

displacement  constraint  case .  4-29 

vii 


4-35.  Reduced  firequency  -  frequency  (k-lB)  curves  for  the  first  4  modes  in  the 

displacement  constraint  case .  4-29 

4-36.  Velocity  -  artificial  damping  (V-g)  curves  for  the  first  4  modes  in  the 

mass  forward  (mg)  case .  4-30 

4-37.  Reduced  frequency  -  frequency  (k-OJ)  curves  for  the  first  4  modes  in  the 

mass  forward  (mg)  case .  4-30 

4-38.  Velocity  -  artificial  damping  (V-g)  curves  for  the  first  4  modes  in  the 

mass  aft  (mg)  case .  4-31 

4- 39.  Reduced  frequency  -  frequency  (k-GJ)  curves  for  the  first  4  modes  in  the 

mass  aft  (mg)  case .  4-31 

5- 1.  Qualitative  view  of  the  effect  of  linking  together  design  variables.  .  5-3 

5-2.  Example  of  possible  displacement  and  natural  frequency  curves  as  a 

function  of  structural  stiffness .  5-7 

a-1.  First  mode  of  vibration  -  stress  constrained  case .  a-2 

a-2.  Second  mode  of  vibration  -  stress  constrained  case .  a-2 

a-3.  Third  mode  of  vibration  -  stress  constrained  case .  a-3 

a-4.  Fourth  mode  of  vibration  -  stress  constrained  case .  a-3 

a-5.  First  mode  of  vibration  -  miiumum  gauge  case .  a-4 

a-6.  Second  mode  of  vibration  -  minimum  gauge  case .  a-4 

a-7.  Third  mod(  of  vibration  -  minimum  gauge  case .  a-5 

a-8.  Fourth  mode  of  vibration  -  minimum  gauge  case .  a-5 

a-5.  First  mode  of  vibration  -  displacement  constrained  case .  a-6 

a-6.  Second  mode  of  vibration  -  displacement  constrained  case .  a-6 

a-7.  Third  mode  of  vibration  -  displacement  constrained  case .  a-7 

a-8.  Fourth  mode  of  vibration  -  displacement  constrained  case .  a-7 

viii 


a-5.  First  mode  of  vibration  -  mass  forward  case .  a-8 

a-6.  Second  mode  of  vibration  -  mass  forward  case .  a-8 

a-7.  Third  mode  of  vibration  -  mass  forward  case .  a-9 

a-8.  Fourth  mode  of  vibration  -  mass  forward  case .  a-9 

a-5.  First  mode  of  vibration  -  mass  aft  case .  a-10 

a-6.  Scoond  mode  of  vibration  -  mass  aft  case .  a-10 

a-7.  Third  mode  of  vibration  -  mass  aft  case .  a-1 1 

a-8.  Fourth  mode  of  vibration  -  mass  aft  case .  a-1 1 


ix 


List  of  Tables 


Table  Page 

3- 1.  Material  Properties . 3-3 

4- 1.  Aerodynamic  Stability  Derivatives .  4-16 

4-2.  Natural  Frequencies  for  Stabilizer  Designs .  4-24 

4- 3.  Flutter  Speed  for  the  6  Design  Cases .  4-26 

5- 1.  Stress  values  <^f  the  element  with  maximum  stress  from  each  of  the  21 

design  variable  groups .  5-4 


X 


A 

[A] 

[B] 
bo 

[C] 
{d} 
{d.) 
{d.} 
{d} 

[D] 

[E] 
{F} 
{F.} 
(F.) 
{F,} 

f. 

[G] 

g 

GW 

H 


List  of  Symbols 

Surface  area  of  Wing  Panel. 

Aerodynamic  Matrix  relating  Force  to  downwash. 

Strain-displacement  matrix. 

Reference  chord  length  (1/2  of  root  chord) 

Damping  matrix. 

Displacement  vector  for  n  discrete  nodal  points. 

Structural  node  displacement  vector. 

Aerodynamic  node  displacement  vector. 

Displacement  modal  vector. 

Downwash-pressure  coefficient  matrix,  or  stress-displacement  matrix. 
Material  property  matrix. 

Force  vector. 

Aerodynamic  force  vector. 

Structural  fence  vector. 

Vector  of  aerodynamic  forces  which  are  not  dependent  on  structural 
displacements. 

Aerodynamic  doublet  line  force  on  panel  i. 

Splining  matrix. 

Artificial  damping  variable. 

Constraint  equation. 

Geometric  properties  of  the  wing,  (size,  sweep,  etc) 

xi 


{IC} 

k 

K 

[K] 

[L] 

li 

[M] 
M 

[N] 
{P} 
P 

P 

q 

[Q] 

{S) 

[S] 

m 

u 

u 

a 


Downwash  across  rigid  wing. 

Reduced  frequency,  (  b,,  CD  /  U  )  or  spring  constant. 

Kernel  Function. 

Finite  element  stiffness  matrix  or  kernel  function  matrix. 

Coefficient  matrix  which  equates  the  value  of  the  incremental 
downwash  due  to  surface  displacements  due  to  flexibility. 

Length  of  bar  element  i. 

Mass  matrix. 

Mach  number. 

Displacement  shape  function. 

Concentrated  forces  at  discrete  points. 

Mcxlal  Pressure  Coefficient. 

Non-dimensional  pressure  coefficient. 

Dynamic  pressure.  (V4  p  U^) 

General  aerodynamic  force  matrix. 

Search  vector  in  design  space. 

Splining  matrix. 

Displacement-downwash  coefficient  matrix. 

Free  stream  velocity 
Displacement  function. 

Velocity  function. 


•  • 
XU 


Acceleration  function. 


Aerodynamic  nodal  point  displacements, 
u,  Structiural  nodal  point  displacements. 

{w}  Vector  containing  the  value  for  downwash  at  n  discrete  points. 

w,w  Downwash  function  across  wing,  modal  downwash  function. 

{X}  Vector  of  design  variables. 

Z  Objective  function. 

C,T1  Integration  constants, 

e  Strain  function  vector. 

0  Stress  function  vector. 

0  Surface  traction  function  or  Lagrangian  function. 

K  Damping  factor, 

p  Mass  density. 

OQ  Frequency  of  oscillation. 

X  Dummy  variable  used  to  calculate  eigenvalues. 

{(|)}  Eigenvector. 


xiii 


AFIT/GAE/ENY/92D-24 


Abstract 


The  purpose  of  this  thesis  is  to  investigate  the  influence  of  structural 
optimization  parameters  on  the  stmctural  and  aeroelastic  behavior  of  a  vertical  tail. 
The  outside  geometry  of  the  tail  is  assumed  fixed.  The  objective  of  the  optimization 
is  to  reduce  the  structural  member  sizes  so  that  the  overall  weight  of  the  tail  is 
minimized.  The  influence  that  design  loads,  design  constraints,  and  design  variable 
linking  has  on  the  optimal  design  is  observed  by  examining  several  design  cases.  The 
static  aeroelastic  properties  and  flutter  speeds  of  the  various  design  cases  are  then 
compared.  The  doublet  lattice  method  is  used  to  calculate  the  steady  and  unsteady 
subsonic  aerodynamic  loads  and  finite  element  theory  is  used  to  model  the  vertical  tail 
structure. 
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AN  INVESTIGATION  INTO  THE  INFLUENCE  OF  STRUCTURAL 
OPTIMIZATION  ON  THE  AEROELASTIC  PROPERTIES  OF 

A  VERTICAL  TAIL 

1.  Problem  Statement  and  Background 

1.1  Problem  Statement 

This  thesis  investigates  the  influence  of  structural  optimization  on  the  the 
minimum  weight  design  and  aeroelastic  response  of  the  support  structure,  or  tailbox, 
of  a  vertical  tail.  The  vertical  tail  is  designed  for  subsonic,  high  angle  of  attack  flight 
conditions.  The  design  loads  are  based  on  the  maximum  aerodynamic  loads  produced 
from  the  tail  during  steady-state  level  flight  The  aerodynamic  loads  and  the  response 
of  the  structure  to  these  loads  are  predicted  with  a  computer  generated  mathematical 
model.  The  exterior  geometry,  or  shape,  of  the  tail  and  tailbox  is  assumed  to  be  fixed. 
Since  structural  weight  is  directly  related  to  the  size  of  the  structural  members,  a 
minimum  weight  design  is  calculated  by  determining  the  minimum  amount  of 
stiuctural  material  which  is  required  to  support  the  design  load  without  violating  any 
of  the  design  constraints. 

To  identify  the  influence  that  design  constraints  have  on  the  structural 
members,  an  optimal  design  is  calculated  for  several  sets  of  design  constraints.  For 
each  constraint  set,  a  minimum  weight  design  is  calculated.  From  the  results  of  each 
of  these  design  cases,  ccmclusions  are  drawn  on  (a)  the  stractural  member  sizes  most 
influenced  by  the  design  constraints,  and  (b)  the  cost  in  structural  weight  for  the 
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desired  aeioelasdc  response. 


1.2  Introduction 

Research  in  the  area  of  aeroelastic  optimization  has  increased  dramatically  in 
the  last  15  years.  This  is  a  result  of  the  emphasis  in  the  engineering  community 
shifting  from  developing  optimal  designs  at  the  test  and  evaluation  level,  back  to 
developing  optimal  designs  at  the  product  design  level.  This  shift  in  interest  has  been 
fueled  by  changes  in  the  aircraft  industry.  Four  of  the  major  changes  are  (35:977) 

(1)  new  aircraft  systems  are  now  very  expensive  to  build  resulting  in 
fewer  aircraft  designs, 

(2)  fewer  current  aircraft  systems  reduce  a  designer’s  ability  to  rely  on 
past  experience  when  making  design  decisions, 

(3)  the  new  aircraft  systems  are  enormously  complex  and  the 
interdependency  of  the  aircraft’s  components  on  the  overall  performance 
is  difficult  for  the  human  mind  to  conceptualize,  and 

(4)  the  processing  capability  of  modem  computers  has  increased 
dramatically  so  as  to  allow  the  designer  to  model  complex  relationships 
numerically. 

These  changes  have  forced  the  designer  to  make  more  and  more  decisions  based  on 
predictions  of  design  behavior  instead  of  actual  test  data.  The  most  common  method 
for  predicting  the  behavior  of  a  design  is  to  create  a  mathematical  model.  These 
models  attempt  to  duplicate  the  characteristics  of  a  proposed  design  and  the  response 
of  the  design  to  given  loading  conditions. 

The  challenge  for  the  engineer  is  to  develop  a  mathematical  model  which  best 
predicts  the  actual  behavior  of  the  design.  If  the  noodel  does  not  closely  duplicate  the 
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behavior  of  an  actual  design,  then  an  optimal  design  based  on  the  mathematical  model 
would  not  be  very  useful.  For  the  case  of  the  vertical  tail,  the  model  usefulness 
depends  on  (a)  the  accuracy  of  the  predicted  aerodynamic  loads  and  (b)  the  accuracy 
of  the  predicted  structural  response  to  the  applied  loads.  Therefore,  before  optimal 
theory  can  be  used  to  design  the  structure  of  the  tail,  the  aerodynamic  and  structural 
models  must  be  examined  to  gauge  the  quality  of  the  optimal  design. 

1.3  Prediction  of  Aerodynamic  Behavior 

Theoretical  prediction  of  aerodynamic  lift  forces  produced  by  lifting  surfaces 
pre-dates  powered  flight  Most  of  the  early  work  focused  on  addressing  the  surface  as 
a  single,  planar  surface  in  steady  flight.  These  classical  approaches  can  be  found  in 
McCormick  (27:44-71)  or  in  most  undergraduate  aerodynamic  texts.  The  prediction  of 
lift  with  these  methods  was  not  very  accurate,  required  empirical  correction  factors, 
and  was  not  applicable  to  oscillating  surfaces.  In  addition,  only  information  on  the 
total  lifting  surface  properties  could  be  calculated.  Information  on  the  local  values  of 
lift,  pressure,  and  velocity  could  not  be  extracted  with  these  methods. 

The  complexity  of  the  governing  equations  restricted  the  bulk  of  research  to 
inviscid,  linear  potential  flow.  In  addition,  the  non-linear  viscous  effects  near  the 
surface  are  ignored  since  their  influence  is  isolated  to  the  surface  and  contribute  little 
to  tail  deflections.  Drag  forces  ate  assumed  not  to  influence  the  aeroelastic  properties 
of  a  lifting  surface  either.  Therefore,  linear  potential  flow  theory  reduces  the 
complexity  of  the  problem  signiflcantly  at  only  a  minor  cost  in  theoretical  accuracy. 
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Theodorsen  published  the  first  major  paper  which  addressed  the  oscillatory 
nature  of  aerodynamics  over  a  lifting  surface  (37:413*440).  Though  this  was  an 
important  step  forward,  his  work  focussed  on  a  two-dimensional  planar  airfoil  which 
oscillated  in  pitch  and  plunge  only  and  did  not  account  for  the  flexible  nature  of  the 
airfoil  or  provide  a  function  of  the  pressure  distribution  across  the  airfoil. 

The  first  significant  success  in  modeling  continuous  pressure  function  across  a 
lifting  surface  is  credited  to  Kussner  (20).  He  developed  a  functional  relationship 
between  the  pressure  on  a  wing  and  the  downwash  at  the  surface  of  the  wing.  The 
relationship  was  developed  from  the  linear  full  potential  equation.  Kussner’s 
relationship  between  pressure  and  downwash  is  referred  to  as  the  kernel  function. 
Unfortunately,  the  kernel  function  is  an  integral  function  which  does  not  have  a  closed 
farm  solatioa  due  to  multiple  order  singularities  in  the  integrand.  Therefore, 
approximations  to  the  kernel  function  have  been  made  in  order  to  obtain  a  working 
relationship  between  downwash  and  pressure. 

A  number  of  researchers  continued  with  Kussner’s  work  by  presenting 
approximations  to  the  kernel  function.  Watkins,  Woolston,  and  Chinningham  (38) 
modified  the  kernel  function  by  first  including  non*planar  effects  and  second 
approximating  terms  in  the  integrand  with  an  exponential  expansion  and  assuming  that 
the  time  varying  pressure  distribution  can  be  represented  by  a  finite  number  of 
pressure  modal  shapes.  This  approximation  presented  the  function  in  a  form  more 
amenable  to  numerical  solution.  Ueda  presented  another  method  of  approximating  the 
integral  with  an  infinite  series  solution  (2:345).  Discrete  approximations  to  the  kernel 
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function  were  used  by  Hedman  for  steady  flow  solutions  to  the  downwash  equation 
(1:280).  Albano  and  Rodden  (1:279-285)  improved  on  Hedman’s  work  by  including 
unsteady  flow  contributions.  Giesing  later  (11:693-702)  included  a  model  which 
would  account  for  interference  panels.  A  thorough  review  of  numerical  lifting  surface 
theory  progress  up  to  1968  was  written  by  Landahl  and  Stark  (21:2049-2061). 

More  recently,  Lottati  and  Nissim  (22:1043-1048)  presented  a  method  of 
integrating  the  kernel  function  in  a  piecewise  fashion  reducing  the  number  of 
calculations  required  for  a  convergent  solution.  Nierkerk  demonstrated  a  method  of 
reducing  the  algebra  and  computation  time  involved  in  solving  the  finite  part  of  the 
integral  (31:1194-1196)  by  using  direct  integration  methods.  Bismarck-Nasr  (2:878- 
880)  exactly  represented  the  kernel  function  as  a  number  of  new  functions.  These  new 
function  have  properties  which  permit  efficient  and  accurate  numerical  evaluation. 

1.4  Prediction  of  Structural  Behavior 

Because  of  the  complexity  of  the  aircraft  structure,  solutions  to  the  structural 
response  based  on  continuum  models  are  not  feasible.  Instead,  a  finite  element 
technique  (6)  is  used.  This  technique  models  the  structure  as  a  number  of  individual 
discrete  parts  or  elements.  The  responses  of  these  elements  to  applied  loads  have 
solutions  which  are  easily  derived.  The  response  of  the  entire  structure  to  an  applied 
load  can  then  be  modeled  as  a  linear  superposition  of  the  response  of  each  individual 
member  or  element  of  the  structure.  This  technique  lends  itself  nicely  to  numeric 
conqrutation.  The  accuracy  of  this  technique  is  governed  by  the  accuracy  of  the 
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element  model  and  the  number  of  elements  used.  Most  finite  element  programs 
available  today  use  elements  which  are  based  on  the  solution  to  the  virtual  work  or 
potential  energy  equation  (6:109-112). 

1.5  Optimal  Theory  Applied  to  Aeroelastic  Design 

Most  of  the  work  in  the  field  of  structural  optimization  as  applied  to  aircraft 
lifting  surfaces  before  the  mid  1970’s  involved  discussion  of  feasible  methods  and 
possible  algorithms  (14:382).  The  algorithms  used  mathematical  programming 
techniques  along  with  the  Hnite  element  method  to  model  a  structure.  These 
algorithms  could  not  be  applied  at  the  time  because  of  the  large  number  of 
computations  which  are  required  to  complete  a  rigorous  optimization.  After  the 
computing  power  of  modem  computers  advanced,  (mid  1970’s)  software  began  to 
emerge.  The  number  of  computations  still  remained  high,  so  researchers  reduced  the 
number  of  calculations  by  relaxing  mathematical  restrictions  or  by  reducing  the 
accuracy  of  the  structural  model.  An  example  of  optimization  theory  applied  to  a 
transport  wing  is  provided  by  Grossman  (12:1050-1058). 

Methods  by  Johnson  (18:438-445)  and  Neill  (30:1021-1027)  use  mathematical 
programming  techniques  to  vary  the  physical  dimensions  of  structural  members  in 
order  to  achieve  an  optimum  design.  This  is  the  simplest  technique  to  apply  since  it 
only  requires  modifications  to  each  of  the  element’s  stiffness  matrices. 
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1.6  Closing 


Much  work  has  been  done  on  developing  the  mathematical  models  which 
predict  the  behavior  of  an  actual  design.  The  accuracy  of  the  aerodynamic  and 
structural  models  are  impressive  as  long  as  the  design  does  not  violate  the  limitations 
imposed  on  each  model.  Since  the  aeroelastic  design  capabilities  are  relatively  new 
only  a  limited  number  of  applications  of  these  models  are  available.  To  date  no 
known  major  system  has  been  designed  using  these  optimal  design  codes.  However, 
As  more  applications  are  published  and  the  results  veritied,  the  trust  in  these  optimal 
codes  should  increase. 

The  results  from  this  thesis  further  builds  on  the  optimal  design  concept 
discussed.  The  variations  in  the  vertical  tail  design  demonstrate  how  optimal  design 
can  be  used  to  influence  aeroelastic  behavior  and  to  identify  which  structural  members 
are  most  influenced  by  the  aeroelastic  constraints. 
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2 


Theory 


2.1  Introduction 

This  chapter  concentrates  on  establishing  the  mathematical  concepts  from 
which  the  results  in  chapter  four  are  based.  These  concepts  include  the  functional 
relationship  between  a  moving  fluid  and  the  pressure  distribution  on  a  lifting  surface, 
and  the  structural  response  to  an  applied  load.  The  aerodynamic  and  stmctural  model 
for  the  tail  are  based  on  these  relationships.  The  characteristics  of  the  flow  field 
around  the  lifting  surface,  structure,  and  applied  loads  are  emphasized.  The  definition 
of  optimization  model  including  the  optimization  goal,  the  optimization  variables,  and 
any  limitations  on  the  variables  is  described.  The  algorithm  used  to  reach  the  optimal 
design  is  also  presented. 

The  aerodynamic  assumptions  of  the  doublet  lattice  paneling  method  (1:279) 
are  given  in  section  2.2.  The  finite  element  stractural  analysis  method  is  discussed  in 
section  2.3.  Section  2.4  develops  the  equations  used  to  tie  together  the  static  or 
dynamic  aeroelastic  loads  to  the  structural  deformations.  The  final  section  addresses 
optimization  theory  as  applied  to  structural  design. 
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2.2  Aerodyiiamic  Theory 


The  aerodynamic  model  used  to  predict  the  loads  on  the  vertical  tail  is  based 
on  the  doublet-lattice  method  presented  by  Albano  and  Rodden  (1:279).  This  method 
assumes  that  the  properties  of  the  flow  field  can  be  modeled  with  linear  potential  flow. 
The  governing  equations  for  linear  potential  flow  are  valid  for  an  inviscid  and 
irrotational  fluid.  In  addition,  it  is  assumed  that  the  fiee-stream  velocity  is  an  order  of 
magnitude  greater  than  any  local  velocity  perturbations  near  the  surface  of  the  vertical 
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Figure  2-1  :  Orientation  of  aerodynamic  velocity  vectors  used  in  equation 
fonnulatitHis. 
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tail.  (i.e.  U«  »  u,,  Uy,  and  u^)  The  free  stream  velocity,  U_,  is  also  assumed  to  be 
time  invariant  in  magnitude  and  direction. 

The  governing  equation  for  linear  potential  flow  is  a  second  order  spatial  and 
time  dependent  equation.  The  form  of  the  equation  includes  the  Prandtl-Glauert 
compressibility  effects.  The  governing  equation  (7:132)  is 


1 

dt^ 

ui- 

dxdt  cfcc^ 


=  0 


(2-1) 


where  V®  =  u.  (u  is  a  vector  of  perturbation  velocities.)  The  equation  must  satisfy 
the  following  boundary  conditions 

as  (2-2) 

^  =  f  ♦  y.f  •  'XwA  a-3) 

Equation  2-3  is  the  surface  boundary  condition.  The  surface  boundary  function  is 
commonly  referred  to  as  the  downwash  function,  w(x,y,z,t).  Note  that  the  downwash 
function  requires  the  flow  to  be  tangential  to  the  surface,  therefore,  the  function  is 
only  dependent  on  the  tail  surface  geometry  and  the  motion  of  the  surface  with  respect 
to  tinoe.  By  multiplying  each  side  of  equation  2-3  by  b^  /  U.  where  b^  is  a  reference 
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(2-8) 


where 

P  =  local  perturbation  pressure.  (Pta»i  -  PJ 
Q.  s  free  stream  density. 
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Both  the  pressure  and  downwash  equations  are  a  function  of  perturbation  potential,  so 
a  functional  relationship  exists  between  downwash  and  pressure.  Since  the  value  of 
the  downwash  function  is  known  for  a  given  tail  geometry,  it  is  desirable  to 
manipulate  equations  2-1  through  2-4,  and  2-8  to  determine  the  differential  surface 
pressme  as  a  function  of  the  surface  downwash.  A  relationship  was  developed  by 
Kussner  (20:14)  by  assuming  that  downwash  and  the  pressure  vary  harmonically  with 
time  or 

MxM)  =  c'"'  ^2-9) 

PCxO'aO  = 

The  functions,  w  and  P,  represent  the  modal  amplitudes  of  the  respective  functions. 
The  resulting  relationship  is  commonly  referred  to  as  the  downwash-pressure  function 
(1:279) 

^  dldr\  (2-10) 

where  k  is  the  reduced  frequency  and  M  is  the  Mach  number.  The  function  K  is  a 
singular  integral  function  referred  to  as  the  kernel  function.  This  function  physically 
relates  the  downwash,  w,  at  the  point  (x,y)  due  to  a  unit  pressure  flux,  P,  at  (^,1)). 

The  limits  of  integration  can  be  reduced  to  the  tail  surface  only  by  the  fact  that  the 
perturbation  pressures  are  defined  to  be  zero  off  of  the  surface.  The  doublet-lattice 
method  is  based  on  this  relationship  between  the  downwash  and  the  pressure  function. 
Due  to  the  singular  nature  of  the  kernel  function,  equation  2-10  caimot  be 
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Figure  2-2  :  Definition  of  kernel  function  coordinate  system. 


inverted  to  supply  the  pressure  distribution  as  a  function  of  downwash.  Therefore,  this 
equation  does  not  have  a  closed  form  solution  for  the  pressure  distribution  across  the 
surface.  Instead,  as  was  discussed  in  the  background  section  of  chapter  one,  the  terms 
in  equation  2-10  must  be  approximated  in  order  to  get  the  equation  in  a  form  which 
provides  a  solution  to  the  pressure  function. 

First,  a  solution  to  equation  2-10  is  evaluated  at  a  number  of  discrete  points 
instead  of  the  continuous  domain.  The  discretization  is  performed  by  modeling  the 
surface  as  a  series  of  planar  trapezoidal  panels.  The  panels  are  arranged  in  columns 
parallel  to  the  firee-stream  velocity  vector.  A  line  distribution  of  acceleration  doublets 


(2:19,130,  and  211)  are  placed  along  the  one-quarter  chord  line  of  each  panel.  The 
doublet  strength  is  assumed  constant  for  each  panel.  The  doublet  strength  is 
determined  by  enforcing  the  known  downwash  conditions  (no  flow  through)  at  the 
mid-span  of  the  three-quarter  chord  line  of  each  panel.  This  point  is  the  downwash 
reference  point.  By  enforcing  the  no-flow  through  condition  at  the  downwash 
reference  point,  reference  1  found  that  the  Kutta  condition  is  empirically  satisfied  at 
the  trailing  edge. 

The  strength  of  the  doublet  lines  of  each  panel  is  represented  by  (1:280) 

(2-11) 

t. 


Doublet  -  Lattice 
Aerodynamic 
Panel  Mesh 


where 


fj  =  constant  vertical  force/unit  length. 

Ij  =  direction  of  panel  one-quarter  chord  span, 
dp  =  differential  length  along  Ij. 

Inserting  this  equation  into  the  down  wash  equation  2-10  reduces  the  surface  integral  to 
a  line  integral  and  simplifies  the  equation  to 


77  =  — fiaxpypX/nUj(ii)]dii 


(2-12) 


(q  =  p  U^)  This  equation  represents  the  downwash  at  a  point  (Xj.yj)  due  to  an 
acceleration  doublet  at  (Xj.yj).  The  total  downwash  produced  at  the  downwash 
reference  point  (Xi,yj)  is  determined  by  calculating  the  influence  of  all  the  panels.  This 
is  simply  the  summation  of  the  contribution  from  each  panel  or 


w  _  A  fj 


h 


(2-13) 


where  n  is  the  number  of  panels. 

Assume  that  the  net  pressure  over  the  panel  is  the  total  force  produced  from 
the  doublets  along  the  panel  span,  divided  over  the  panel  area  or 


?  =  =  _ 4 _  (2-14) 

Sj  AXy  /y  COSAy  AXy  COSAy 

where 

Ax,  =  average  panel  chord  length. 

Aj  s  panel  one  quarter  chord  span  sweep  angle. 

By  substituting  these  equations  into  the  downwash  equation  2-10,  the  downwash 
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equation  simplifies  to 


(2-15) 


where  q  =  Vi  and  the  integral  term  is  evaluated  numerically.  A  vector  of 
downwash  values  can  then  be  calculated  and  the  resulting  equations  can  be  written  in 
matrix  form  * 


inverting  [D] 


< 


=  [D]-' 


(2-16) 


(2-17) 


Equation  2-17  then  provides  the  pressure  at  each  of  the  n  panels  given  the  downwash 
values  at  the  n  points.  The  overall  force  and  moment  values  are  determined  with  an 
integration  matrix  [S]  (23:86).  This  matrix  acts  on  [D]'*  (wAJ.)  and  transforms  tiie 
pressures  into  forces  and  moments.  Before  integrating  the  pressures,  each  side  of 
equation  2-17  is  multiplied  by  q  (dynamic  pressure)  so  that  the  pressure  values  are 


Throughout  this  thesis  a  single  solid  line  inside  of  vector  notation,{— },  denotes 
division,  whereas  a  double  line  (-  -}  denotes  vector  partitioning. 
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dimensional.  The  total  lift  and  moment  about  a  given  reference  point  then  is 
represented  with 


F 

M 


= 


-1 


w 

U. 


(2-18) 


The  aerodynamic  forces  and  moments  produced  by  the  tail  can  now  be 
calculated  for  a  given  tail  geometry.  The  fmces  are  based  on  linear  potential  flow 
which  varies  harmonically  with  respect  to  time  across  the  tail.  Steady  flow  results  are 
obtained  by  setting  the  frequency,  C3,  equal  to  zero.  The  linear  flow  assumption 
restricts  the  analysis  to  small  angles  of  attack  because  the  perturbation  velocities  must 
remain  much  less  than  U,..  Finally,  it  is  assumed  that  the  pressure  variation  can  be 
modeled  with  a  finite  number  of  acceleration  doublet  lines  instead  of  a  continuous 
sheet  of  doublets.  How  these  forces  are  used  in  conjunction  with  the  structural  model 
will  be  explained  in  the  aeroelastic  model  section  (2-4). 


2.3  Structural  Analysis 

The  finite  element  method  (FEM)  is  used  to  model  the  structural  response  of 
the  tail  to  an  applied  load.  This  method  is  a  numoical  procedure  for  analyzing 
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complex  structural  behavior.  Instead  of  modeling  an  entire  structure  with  the 
governing  continuum  equations,  the  structure  is  modeled  as  an  assemblage  of  discrete 
pieces.  These  pieces  are  referred  to  as  elements.  The  behavior  of  each  of  the 
elements  is  relatively  simple  to  model  with  available  equations  of  motion,  (i.e.  beam 
bending,  torsion  bar,  plane  stress,  etc.)  llie  method  is  based  on  satisfying  the 
governing  equations  at  a  discrete  number  of  points  and  assuming  the  form  of  the 
displacement  between  these  points  (interpolation).  The  elements  are  defined  by 
specifying  the  connectivity  of  the  element  to  these  points,  called  nodal  points.  Each 
nodal  point  has  up  to  six  degrees  of  motion  (3  translations,  3  rotations).  The 
individual  elements  are  intercoimected  through  the  nodal  points.  The  compatibility 
and  equilibrium  constraints  on  the  governing  equations  are  imposed  at  these  points. 

To  simplify  the  governing  equations  the  following  assumptions  are  made 

1.  strain  (e)  is  assumed  to  be  linear  widi  respect  to  stress  (a). 

2.  The  structure  deforms  elastically,  (no  residual  strain). 

3.  Damping  is  neglected  (conservative  system). 

4.  Applied  forces  are  not  a  function  of  displacement. 

The  govenung  equation  for  a  finite  element  model  with  damping  is  (6:75) 


[MI  [Jt  ♦  [A]  [di  ♦  m{d)=  ino) 


(2-19) 


where 

[M]  s  Model  mass  matrix. 

[B]  =  Model  damping  matrix. 

[K]  =  Model  stiffness  matrix. 

(d)  s  Generalized  nodal  displacement  vector, 

(dj  =  Generalized  nodal  velocity  vector. 

(d)  -  Generalized  nodal  acceleration  vector. 

{F(t))  s  Nodal  face  vector. 
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which  is  derived  from  the  virtual  work  equation.  Since  damping  is  neglected,  the 
equation  can  be  simpliried  to 

[M]  +  [Kl  id]  =  {fXr)}  (2-20) 

The  displacement  at  any  point  (x,y,z)  in  the  structure  is  based  on  the  respective 
displacement  at  the  nodal  points.  In  other  words 

{«(W))  =  [N1  id]  (2-21) 

where  {u}  represents  the  continuous  displacement  Held  and  [N]  is  the  shape  matrix 
which  acts  as  the  interpolation  matrix  with  the  individual  nodal  points. 

If  the  inertial  contribution  ( [M](d}  )  is  small,  equation  2-20  simplifies  to 

ra  {d\  -  (F)  <2-22) 

The  displacement  at  each  of  the  nodal  points  can  then  be  determined  by  inverting  the 
stiffness  matrix 

{<0  =  [«]-*  (F)  (2-23) 

Therefore,  given  a  set  of  applied  loads  and  the  stiffness  properties  of  the  finite  element 
model,  the  resulting  displacements  can  be  calculated. 

The  stress  and  strain  can  be  calculated  from  the  known  displacements  by  using 
the  strain-displacement  relationship 

le)  -  -  lemidi  -  mw  <2-24) 

where  [d]  is  the  derivative  matrix  which  transforms  displacements  to  strain,  and  [B]  is 
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referred  to  as  the  strain-displacement  (6:17)  matrix  ( [B]  =  [9]  [N]  ).  The  stress-strain 
relationship  is  {o}  =  [E]  {e}  where  [E]  is  the  material  stiffness  properties  matrix 
(6:18).  By  incotporating  the  results  from  equation  2-24 

{0}  =  mimd)  (2-25) 

If  the  inertial  effects  can  not  be  neglected,  then  the  mass  term  and  nodal 
acceleration  term  must  be  included  into  the  solution.  To  simplify  the  resulting  2nd 
order  differential  equation,  the  time  dependency  of  the  displacement  vector  is  assumed 
to  be  harmonic 

{d\  =  {d]  e*"'  (2-26) 

After  substituting  this  into  equation  2-20 

(-«*  [Ml  +  =  {fXOl 

The  model  response  becomes  infinite  for  values  of  o)  where  the  determinant  of 
I  -a>*  [M]  +  [K]  I  is  zero.  This  type  of  problem  is  referred  to  as  an  eigenvalue 
problem.  The  solution  involves  calculating  values  of  where  the  determinant 
vanishes.  These  values  for  oD)  can  be  substituted  back  into  the  governing  equation 
({F)  =  0)  to  obtain  the  corresponding  modal  shape  vectors, 

(^Um]-[K\)  =  0.  (2-28) 

The  mode  shapes  are  referred  to  as  eigenvectors.  These  eigenvectors  are  stored  in  a 
square  matrix  of  n  eigenvectors, 

Since  the  model  is  linear,  the  response  to  any  set  of  applied  loads  is  simply  a 


2-13 


linear  superposition  of  each  of  the  n  modal  responses  fi(t).  For  free  vibration,  the 
response  is 

{d(0)  =  Efl,  {♦,)  (2-29) 

where  the  coefficient  a,  is  the  amplitude  of  the  ith  mode. 

The  preceding  paragraphs  demonstrated  the  usefulness  of  the  FEM  in 
determining  the  structural  response  of  the  vertical  tail  to  both  static  and  time 
dependent  loads.  To  arrive  at  equation  2-23  and  2-27,  the  assumptions  had  to  be 
made  that  the  structure  is  continuous  in  a  piecewise  manner,  the  force-displacement 
relationship  is  linear,  and  the  material  deforms  elastically.  The  forces  can  only  be 
applied  at  the  structural  nodal  points  so  an  interpolation  scheme  is  employed  to 
transfer  the  loads  calculated  in  section  (2.2)  to  the  structural  nodal  points.  The  method 
is  explained  in  the  following  section.  Also,  the  aerodynamic  forces  are  a  function  of 
displacement  which  requires  special  attention.  This,  too,  is  explained  in  the  following 
section. 


2.4  Aeroelastic  Analysis 

This  section  develops  the  equations  which  tie  together  the  aerodynamic  and 
structural  model.  Two  different  approaches  are  developed;  the  first  is  for  the  static 
aerodynamic  case,  and  the  second  is  for  the  critical  dynamic  response,  or  flutter,  case. 
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The  important  issue  with  aerodynamic  loads  is  that  the  loads  are  dependent  on  the 
local  deflections.  The  equations  in  section  2.3  assumed  that  the  loads  were  not 
dependent  on  displacements.  The  load  dependency  on  displacement  results  in  an 
aerodynamic  ’stiftiiess’  term  which  is  incorporated  into  the  structural  stiffness  term. 

Static  Aeroelasticity:  For  the  static  aeroelastic  case,  the  structural  force  equilibrium 
equation  given  in  section  2-3  is 

IK\  [d]  =  {F]  (2-30) 

Assume  that  the  force  vector  is  solely  produced  from  aerodynamic  loads  and  can  be 
divided  into  two  force  vectors 

{FMF,}  *  [F^]  (2-31) 

where  {F,}  is  a  vector  of  aerodynamic  forces  which  are  a  function  of  the  structural 
displacements  and  {Fj}  is  a  vector  of  forces  produced  from  the  tail  as  if  it  were  rigid. 
The  vector  {F,}  is  assumed  to  be  zero  for  zero  nodal  displacements  (rigid)  and  the 
value  for  {Fj}  is  determined  from  the  initial  conditions  (tail  angle  of  attack).  Assume 
the  dependence  of  (Ft)  on  the  structural  displacements  can  be  represented  by 

=  q  [A]  [d]  (2-32) 

where  q  is  dynamic  pressure  and  [A]  acts  as  an  aerodynamic  stiffness  quantity. 
Equation  2-30  can  then  be  rewritten 
or 
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m  -  q[A])  [d]  =  {Fj} 


(2-33) 


[K^{d)  =  {Fj}  (2-34) 

Remember  that  (Fj)  is  independent  of  the  deflections.  If,  as  q  is  increased,  the 
determinant  of  [KT  decreases  in  value  towards  zero;  the  displacements,  {d},  will  have 
to  increase  to  retain  the  equality  of  equation  2-34.  If  the  increase  in  displacement  is 
only  minimal,  then  there  is  no  need  for  concern,  but  if  the  increase  in  nodal 
deflections  is  large  then  the  structure  could  exceed  elastic  limits.  This  phenomena  is 
called  divergence. 

There  is  one  major  problem  with  equation  2-31  through  2-34,  the  equations 
require  that  the  forces  be  applied  at  the  structural  nodal  points,  {d,}.  In  the  case  of 
this  aerodynamic  model,  the  forces  are  located  at  aerodynamic  nodal  points,  (d,}.  The 
aerodynamic  loads  can  only  be  utilized  in  the  structural  model  if  they  act  at  the 
structural  nodal  points.  Therefore  an  interpolation  scheme  called  ’splining’  (23:82)  is 
used  to  transform  the  aerodynamic  frame  of  reference  into  the  structural  frame  of 
reference.  First,  a  matrix,  [G^,  is  developed  which  maps  the  structural  deflections 
into  the  aerodynamic  frame  of  reference.  In  other  words  [G,_J  calculates  the 
equivalent  aerodynamic  displacements  from  the  given  structural  displacements.  The 
interpolation  equation  (23:82)  is 

+  fljX  +  Ojy  +  P^  (2-35) 

where  K,(x,y)  =  (l/ltin)  r*  In  r* 
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CASEl 


■  Nodal  Points  of  known 
values 


CASE  2 


Figure  2-4  :  Two  1-D  examples  of  inteipoladon  and  extrapolation  curves  determined 
from  equation  2-35,  given  three  data  points. 

r,*  =  (x  -  xf  +  (y  -  yj)* 

aj ,  P,  are  coefficients  which  satisfy 
I  Pi  =  I  Xj  Pi  =  Z  yi  P,  =  0. 

The  coefficients  are  determined  by  repeating  this  equation  for  the  displacements  at  the 
n  structural  nodal  points,  (xj.yj),  to  form  a  matrix  equation 


(2-36) 


where  [Q]  contains  the  constants  for  the  structural  nodal  points.  The  [CJ  matrix  is 
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then  inverted  to  find  the  coefficients.  The  displacements  at  the  aerodynamic  grid 
points  are  then  represented  by 


{dj  =  [CJ  [Cy 


(2-37) 


where  [CJ  contains  the  constants  for  the  aerodynamic  nodal  points.  Interpolation  and 
extrapolation  is  then  done  for  the  aerodynamic  displacements.  These  displacements 
will  follow  the  trend  of  the  structural  nodal  points.  If  the  structural  nodal  points 
deflect  in  a  parabolic  pattern  extrapolation  will  continue  the  pattern.  If  the  structural 
nodal  points  deflect  linearly  (for  example,  wing  twist),  the  leading  and  trailing  edge 
aerodynamic  panels  will  deflect  in  the  same  linear  fashion.  As  the  extrapolation  is 
used  for  data  points  further  from  the  stmctural  nodal  points,  the  aerodynamic  point 
deflections  will  tend  to  converge  to  a  linear  behavior. 


In  simpler  form,  equation  M^tet|  ^ 


(2-38) 


To  insure  that  the  overall  force  applied  to  the  structural  model  is  identical  to  the  force 
applied  to  the  aerodynamic  model,  identical  virtual  work  of  both  systems  is  enforced 
or 


(M.)’'  W'.l  -  IH)''  I'".) 

where  the  symbol  5  refers  to  the  virtual  values  of  displacement,  and  (F.)  and  (F,)  are 
the  forces  in  the  aerodynamic  frame  of  reference  and  the  structural  frame  of  reference, 
respectively.  Equation  2-38  and  2-39  can  be  combined  to  obtain 
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Since  virtual  displacements  are  arbitrary  then  the  terms  inside  the  bracket  must  equal 


zero  and 
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This  can  now  be  used  in  equation  2-18  which  results  in 
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(2-42) 


and  provides  the  aerodynamic  loads  at  the  structural  nodal  points. 

In  order  to  determine  the  contents  of  the  matrix  [A]  from  equation  2-32, 
downwash  must  be  presented  as  a  function  of  the  displacement  vector.  This  is  done 
with  a  derivative  matrix  [T]  so  that 


or 


(2-43) 


■^1  = 


(2-44) 


Equation  2-7  illuminates  the  task  of  this  matrix.  The  [A]  matrix  then  has  the  form 
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[A]  =  «[G,J^[S][D]-'[71[G,J  (2-45) 

The  only  unknown  now  left  in  equation  2-34  is  the  structural  displacement 
vector  {d}.  A  solution  to  the  structural  displacements  are  then  obtained  by  inverting 
the  modified  stiffness  matrix  [K^.  These  displacements  can  then  be  used  to  calculate 
member  stresses,  strains,  etc. 

Flutter  analysis:  Aerodynamic  forces  naturally  will  vary  due  to  wind  gusts  or 

maneuvers.  The  wing  responds  to  the  fluctuation  of  aerodynamic  loads  by  oscillating. 
At  low  dynamic  pressures,  the  oscillatory  aerodynamic  loads  typically  act  opposite  to 
the  motion  of  the  wing.  This  effectively  dampens  the  response  of  the  wing.  As  the 
dynamic  pressure  increases,  the  oscillatory  aerodynamic  forces  may  begin  to  dampen 
the  oscillations  less  and,  in  fact,  begin  to  reinforce  the  motion  of  the  wing.  This 
reinforcement  can  quickly  cause  wing  components  to  exceed  elastic  limits  and  fail. 

The  motion  of  the  wing  is  modeled  using  equation  2-17  of  section  2.2.  As  in 
divergence,  the  main  difficulty  is  modeling  the  aerodynamic  forces  accurately. 
Assuming  no  structural  damping  and  concentrating  on  the  oscillatory  aerodynamic 
forces  only,  the  governing  equation  for  the  structural  response  can  be  written 

=  qKXMXm  [d] 

where  [Q]  is  the  aerodynamic  stiffness  matrix  and  the  other  terms  were  defmed 
previously.  With  oscillatory  motion,  [Q]  is  a  complex  matrix  which  introduces 
danqnng  into  the  problem.(  the  dancing  terms  come  from  the  kernel  function) 
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Before,  a  solution  to  this  second  order  differential  equation  involved  assuming  the 
motion  was  harmonic  and  oscillated  at  a  frequency  e‘“  To  model  this  problem  as 
harmonic  then  the  damping  from  [Q]  must  be  removed.  To  artificially  remove  the 
damping  from  the  problem,  an  artificial  damping  factor  (g)  is  inserted  into  the 
equation.  If  a  value  for  g  is  chosen  which  cancels  out  the  aerodynamic  damping  of 
one  of  the  modes  then  the  system  can  oscillate  harmonically  in  that  mode.  The  sign 
of  the  damping  factor  is  important.  When  the  sign  of  g  is  negative,  the  aerodynamic 
damping  is  stable  and  would  have  caused  the  oscillations  to  die  out.  If  g  is  positive, 
then  the  aerodynamic  damping  would  have  caused  unbounded  growth  of  the  response. 
A  velocity/Mach  number  combination  at  which  g  is  zero  reflects  a  point  at  which  there 
is  no  aerodynamic  damping  in  one  of  the  modes.  This  is  the  critical  point  at  which 
one  of  the  modes  of  vibration  become  unstable.  The  smallest  velocity/Mach  number 
pair  for  which  g=0  is  considered  to  be  the  flutter  point. 

Equation  2-46  is  rewritten  to  reflect  the  harmonic  assumption  and  to  include 
the  artificial  damping  term 

(  -[A<]a)'+(l+<g)[A]-9[<21){rf,)«‘"*=0.  (2-47) 

where  {d,}  is  the  modal  participation  vector.  The  solution  to  this  equation  is  solved 
with  a  technique  referred  to  as  the  K  method  (23:2.6-1).  Equation  (2-47)  is  rewritten 
to  accommodate  this  method  and  the  steps  follow.  Combine  [Q]  and  [M]  within  one 
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set  of  parenthesis 


[<a 

m 

2{2k) 

l+ig 

WJ=0. 


(2-48) 


The  equation  is  rewritten  by  factoring  out  the  (c/2k)^  term. 


{dj=0. 


With  this  representation,  the  eigenvalue  is 


-if 

(a+,P)2=^ 

1+^g 


or 


(2-49) 


(2-50) 


(2-51) 


(«+iP)  =  +  iU^ 

From  this,  the  velocity  and  damping  factor  can  be  found  to  be 


«.=P 


(2-52) 


(2-53) 


Note  that  as  g  0,  the  imaginary  part  of  eigenvalue  also  tends  to  zero.  Solution  to 
equation  2-49  dien  requires  the  user  to  input  a  value  fn*  k,  m,  and  p.  Inputting  these 
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values  into  the  equation  and  solving  the  resulting  eigenvalue  problem  will  produce  a 
solution  for  g  and  V.  The  flutter  speed  is  the  velocity  for  which  the  artificial  damping 
factor  is  zero.  For  a  more  complete  investigation  on  this  subject,  the  reader  is  referred 
to  references  3  and  5. 

2.5  Structural  Optimization 

The  optimal  design  for  the  vertical  tail  is  determined  using  a  numerical 
optimization  algorithm.  This  section  introduces  the  general  problem  statement  and  the 
basic  concepts  and  definitions  of  the  optimization  problem.  The  application  of  these 
concepts  to  the  structural  design  problem  is  then  outlined. 

Before  an  optimal  design  can  be  determined  three  characteristics  need  to  be 
known  about  the  design;  (1)  what  end  result  is  desired,  (2)  what  changes  in  the  design 
are  allowed  to  best  attain  that  result,  and  (3)  what  limitations  are  there,  to  these 
changes.  To  mathematically  model  the  design  process  each  of  these  characteristics  is 
rqmsented  as  function  relationships.  The  first  charactraistic  is  called  the  objective  of 
the  design.  The  values  that  are  allowed  to  be  changed  during  the  optmiization  process 
are  referred  to  as  the  design  variables.  The  objective  function  should  be  a  function  of 
the  design  variables.  This  way  the  design  variables  can  be  changed  in  order  to  better 
meet  the  design  objective.  The  objective  function  is  constructed  so  that  the  optimal 
design  is  reached  only  when  die  function  is  at  it’s  minimum  value.  The  design 
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variables  most  often  are  restricted  to  a  range  of  values.  If  the  value  of  all  of  the 
design  variables  is  within  this  defined  range  then  the  design  is  considered  to  be  a 
feasible  design.  If  any  of  the  design  variables  fall  outside  of  the  defined  range  the 
design  is  considered  to  be  an  infeasible  design.  The  range  is  defined  by  an  equation 
called  the  constraint  equation.  A  design  may  have  many  constraint  functions. 

For  the  stmctural  design  problem  the  objective  function  could  be  to  minimize 
the  structural  weight  by  decreasing  the  sizes  of  all  the  structural  members.  If  the 
weight  was  directly  proportional  to  the  member  sizes  and  the  problem  defiiution  did 
not  state  any  constraints,  it  should  be  obvious  to  see  that  the  optimal  design  is 
achieved  when  the  member  sizes  are  zero.  This,  of  course  is  not  true  so  constraints 
are  placed  on  the  member  sizes  to  account  for  maximum  stress  values,  maximum 
displacement  allowed,  and  any  other  structural  behavior  that  the  designer  wishes  to 
avoid.  The  finite  element  method  is  used  to  produce  the  equations  necessary  for  the 
constrained  optimization  problem.  The  optimal  solution  is  determined  with  an 
iterative  scheme  referred  to  as  the  method  of  modified  feasible  directions  (20: A- 11). 

The  objective  and  constraint  function  can  be  non-linear,  but  each  function  and 
its  first  derivative  must  be  continuous  and  be  defined  for  every  point  in  the  design 
space.  The  general  non-linear  constrained  optimization  problem  is  stated  as  follows; 

Find{X)=LX,.Xa,....X.J 

such  that 

Z  =  F({X})  Minimum  Value, 
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for  values  of  {X}  which  satisfy  the  constraints 
Gj(X)  ^0  j  =  1,2......C 

Xi'“*~  ^  Xi  <  i  =  i,2,....,n 

where  F({X})  is  the  objective  function,  {X}  is  a  vector  of  the  design  variables,  and 
G({X})  is  the  constraint  function.  As  an  example  of  a  two  design  variable,  (X„  Xj), 
constraint  problem,  consider  a  constraint  function  G({X})  s  a^  -  (Xj^  +  Xi^).  (see 
figure  4)  A  special  case  of  the  constraint  equation  is  the  side  constraint.  This 
constraint  limits  directly  the  value  of  a  design  variable,  such  as  X,  >  0  and  Xj  >  0. 
(see  figure  4)  Before  explaiiung  how  these  equations  can  be  applied  to  structural 
optimization,  an  overview  of  how  the  solution  to  this  problem  is  found  is  necessary. 

The  solution  of  this  problem  is  determined  with  an  iterative  first  order 
procedure  called  the  method  of  modified  feasible  directions.  An  initial  set  of  design 
variable  values,  (X"},  must  first  be  given.  Assume  that  these  values  do  not  violate 
any  of  the  constraints.  This  is  referred  to  as  a  feasible  design.  A  new  guess  on  an 
improved  design  is  then  calculated  by 

+  a  (2-54) 

where  q  is  the  iteration  number,  {S**}  is  a  vector  which  defines  the  direction  in  the 
design  space  where  the  new  guess  will  rrx>ve,  and  a  defines  the  magnitude  or  distance 
of  the  move.  Values  for  the  search  vector,  (S),  and  a  must  now  be  determined. 

The  search  vector  is  required  to  point  in  a  direction  which  reduces  the  value  of 
the  objective  function,  F({X)), 
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Figure  2>5:  Example  optimization  design  space.  (F(X)  objective  function) 


VF  •  {5}  s  0 


(2-55) 


or 


I  VF  I  •  \{S]\  cose  i  0  (2-56) 

The  gradient  of  F({X})  points  represents  the  direction  for  which  F  increases  most 
rapidly,  therefore,  the  direction  for  most  rapid  decrease  is  the  negative  of  that 
direction.  As  long  as  the  cosine  of  the  angle  between  the  two  directions  is  less  than 
zero,  then  {S}  is  in  a  direction  which  will  decrease  the  objective  function. 

For  the  starting  vector,  (X”),  assume  that  the  values  are  not  close  to  any 
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constraint  boundaries.  For  this  case,  the  direction  which  would  minimize  the  objective 
function  most  quickly  is  the  negative  of  the  gradient  of  F({X}),  or  -VF.  So  for  the 
unconstrained  problem  the  search  direction  would  be 


[S]  =  - 


VF 

|VF1 


(2-57) 


where  {S}  is  a  normalized  with  respect  to  the  euclidean  norm  of  VF.  Now,  the 
magnitude  or  distance  of  the  search  vector,  a,  must  be  determined. 

The  value  of  a  is  limited  to  a  magnitude  which  does  not  cause  one  of  the 
constraint  equations  to  be  violated.  To  determine  a  value  for  a,  the  constraint 
functions  are  calculated  at  a  number  of  prescribed  values  for  oc.  When  the  first  value 
of  a  for  which  a  constraint  is  violated,  Gj({X})  >  0,  is  calculated,  an  polynomial 
interpolation  scheme  is  used  to  estimate  the  maximum  value  of  a  which  can  be  used 
which  will  not  violate  the  constraint.  This  value  is  then  used  for  the  iteration 
equation. 

By  using  equation  2-54,  a  new  guess  value  for  the  design  variables  can  now  be 
calculated.  This  value  should  be  close  to  one  or  more  of  the  constraint  boundaries. 
The  new  vector  (X)  will  result  in  a  reduced  value  of  F((X}). 

The  new  vector  {X**}  still  might  be  able  to  be  changed  so  that  the  objective 
function  is  reduced  even  though  it  is  against  one  or  more  constraint  botmdaries.  First, 
the  new  search  direction  must  still  satisfy  equation  2-55  but  to  insure  the  direction 
does  not  violate  the  constraint  equations,  the  search  direction  must  also  satisfy  the 
constraint 
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VG>({X})  •  {S}  i  0 


(2-58) 


which,  as  in  equation  2-56,  also  states  that  the  angular  difference  between  the  VG  and 
{S}  must  be  between  90  and  270  degrees. 

If  during  the  optimization  process,  the  design  variable  vector  moves  to  a  point 
in  the  design  space  that  one  or  more  constraints  are  violated,  equation  2-55  and  2-58 
are  modified  slightly  to  bias  the  (S)  vector  in  a  direction  which  will  return  {X}  to  the 
feasible  region  (20:77). 

The  process  is  considered  to  have  converged  when  the  objective  function  can 
not  be  reduced  further.  For  example,  when  a  modified  design  does  not  reduce  the 
objective  function  by  0.1  percent.  Several  other  convergence  criteria  are  used  to 
insure  that  the  code  does  not  continue  unnecessarily.  These  criteria  are  useful  in 
problem  specific  cases,  so  they  are  not  discussed  here. 

This  was  a  quick  overview  of  the  method  of  modified  feasible  directions.  For 
the  sake  of  brevity  the  mathematical  intricacies  as  well  as  the  variations  of  this  method 
for  particular  cases  has  been  excluded.  However,  a  thorough  explanation  can  be  found 
in  references  14,  23,  or  28.  The  attention  will  now  be  shifted  to  the  structural  design 
problem  and  how  the  equations  of  force  equilibrium  are  modified  in  order  to  be  in  a 
form  similar  to  the  constrained  optimization  problem. 

The  design  properties,  or  design  variables,  for  the  vertical  tail  are  the  cross- 
sectional  area  and  thickness  of  the  wing  box  bars  and  panels  respectively.  It  is 
assumed  that  the  matoial  is  homogenous  so  that  the  weight  varies  directly  with  the 
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Figure  2-6:  Graphical  representation  of  VF,  VG,  and  (S). 

material  volume.  The  value  for  these  properties  must  then  be  reduced  in  order  to 
produce  a  lighter  weight  design.  The  design  variables  effect  the  force  equilibrium 
equation  through  the  stiffness  matrix  [K]  and  the  mass  matrix  [M].  For  example  the 
stiffness  matrix  for  a  simple  two  degree  of  freedom  rod  has  a  common  factor  of  EA/L 
where  E,  A,  and  L  represent  the  modulus  of  elasticity,  cross-sectional  area,  and  bar 
length  respectively.  So,  the  stiffness  matrix  can  be  written  as  [K]  =  EA/L  FK^^I 
(23:125)  where  A  represents  the  design  variable  and  is  the  factored  [K]  matrix. 
For  the  mass  matrix  of  the  bar  element,  the  cross-sectional  area  can  be  factored  out  of 
the  mass  matrix  so  that  [M]  =  p  A  L  [MiaoJ  where  p  is  the  matoial  mass  density  and 
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[Hnodl  is  the  factored  mass  matrix.  Similar  operations  can  be  done  on  other  element 
stiffness  and  mass  matrices  although  the  design  variable  may  not  be  factorable  out  of 
the  respective  matrices.  The  design  variables  now  are  used  to  set  up  the  objective  and 
constraint  function. 

The  objective  of  the  design  is  to  minimize  the  weight  of  the  overall  structure. 
The  weight  of  the  design  can  be  determined  by  summing  the  weight  of  each  of  the  m 
elements  or 


F({XJ)  =  Weight  =  f:  (2-59) 

i 

where  V(Xi)  is  the  volume  of  the  ith  element  For  a  bar  element  the  volume  is  V(Xi) 
=  jc  AjL/4.  The  variable  Lj  is  fixed  since  geometry  of  structure  is  given. 

The  form  of  the  objective  function  for  the  stmctural  design  problem  is  now 
known.  It  is  desirable  to  modify  the  design  variable  definition  slightly  to  linearize 
some  of  the  constraint  equations.  The  design  variable  is  defined  as 

X,  =  -1  (2-60) 

The  advantage  gained  by  this  definition  will  be  seen  during  the  constraint  equation 
develt^ment. 

Assuming  the  design  variables  are  the  cross-sectional  area,  the  gradient  of  the 
objective  function  in  equation  2-59  is 

This  can  then  be  used  to  calculate  the  direction  of  steepest  descent 

Several  constraints  are  used  during  the  design  of  the  vertical  tail.  In  particular. 
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(2-61) 


stress  constraints,  displacement  constraints,  frequency  constraints,  and  size  constraints 
are  used  to  limit  the  values  of  feasible  design  variable  values.  The  constraint 
equations  as  well  as  the  constraint  gradient  which  is  used  in  equation  2-58  are 
presented  next. 

A  minimum  size  constraint  is  a  linear  constraint  which  is  also  referred  to  as  a 
side  constraint.  It’s  purpose  is  to  limit  the  upper  and/or  lower  bound  of  the  design 
variable.  This  constraint  might  be  enforced  because  of  manufacturing  limitations  or 
assembly  and  fit  restrictions.  The  constraint  equation  is 

W  =  ^  0  2-62) 

=  iX^-X)  s  0 

Where  and  X,^  are  the  respective  limits  of  the  design  variables.  There  would 
be  2  equations  for  each  of  the  m  design  variables  and  the  gradient  of  each  of  these 
equations  would  be  a  vector  of  zeros  except  at  the  ith  entry  which  would  have  a  ±1 
depending  on  whether  it  was  the  upper  or  lower  limit  equation. 

For  a  displacement  sensitive  design,  a  displacement  constraint  equation  would 
restrain  the  displacement  at  a  node  j  to  a  minimum  and  or  maximum  value.  The 
constraint  equation  is  similar  to  the  minimum  gauge  function  and  is  defined  as 
where  the  subscript  j  denotes  the  jth  node.  Note  that  the  constraint  equation  will  now 
most  likely  be  a  function  of  most  or  all  of  the  design  variables.  The  gradient  function 
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GiW)  (2-6^ 

G({X])  =  (d^  -  dp^O 

for  displacement  constraints  represent  the  change  in  displacement  due  to  a  change  in 
each  of  the  design  variables. 


VG= 


-Jt 

am 


lower  limit 
upper  limit 


(2-64) 


The  gradient  function  is  derived  from  the  force  equilibrium  equation  (29:134).  The 
force  equilibrium  equation  developed  in  section  2.3  is 


lK\{d}={F]  (2-65) 

Let  the  design  variable  for  each  element  be  represented  by  {X}.  The  change  in 
equation  2-58  due  to  a  change  in  the  design  variable  can  be  represented  by 


Assume  that  (F)  does  not  vary  with  respect  to  Xj.  Then  equation  2-45  can  be  written 
as 


(2-67) 

All  the  variables  on  the  right  hand  side  of  the  equation  are  known  from  the  problem 
formulation  and  initial  analysis.  This  equation  provides  the  gradient  of  the 
displacement  constraint  function  with  respect  to  any  of  the  design  variables  (X). 
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The  stress  constraint  equation  limits  the  values  of  the  design  variables  so  that 
the  yield  stress  values  for  each  element  are  not  violated.  Therefore  there  will  be  a 
constraint  equation  for  each  element.  The  constraint  equation  is  written 

GfiX)  =  ^  0  (2-68) 

~  ^taax-ComprtMiioH  ^  0 

where  a„„  is  the  tensile  or  compressive  yield  stress.  The  governing  equation  for 
member  stress  was  given  by  equation  2-25  and  is  repeated 

{o\=mm{d}  (2-69) 

The  dependence  of  stress  on  the  design  variables  is  determined  by  taking  the 
derivative  of  {a}  with  respect  to  Xj  which  is 

■>  P-W) 

where  [D]  =  [E]  [B] 

Substituting  in  equation  2-67  and  noting  that  =  [0.0] 

oX^ 

(2-71) 

0A|  CMkj 

where  all  the  tnins  on  the  right  hand  side  are  known  from  the  initial  problem  analysis. 

The  last  constraint  that  is  enforced  on  the  tail  design  is  a  frequency  constraint 
This  constraint  requires  that  the  first  natural  frequency  be  less  than  some  given 
maximum  or 
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G({X})  =  a>f-o>Lx  ^  0 


(2-72) 


This  constraint  would  put  a  ceiling  on  the  maximum  value  for  the  fundamental 
frequency.  This  might  be  important  to  the  functioning  of  electronic  equipment  or  the 
fatigue  life  of  a  component.  Note  that  is  used  instead  of  co  to  simplify  the  gradient 
function.  The  gradient  of  this  function  then  is 


VG({X})  = 


dcij 

aw 


(2-73) 


The  gradient  of  the  frequency  is  determined  from  equation  2-28  by  noting  that  for  the 
ith  frequency 


-  «f[Ml{4),l=0  (2-74) 

where  4>,  is  the  ith  mode  shape  or  eigenvector.  Differentiating  this  equation  with 
respect  to  (X,}  results  with 


2  dcji 


(Mil*,)  -  - 


-O, 


(2-75) 


If  it  is  assumed  that  the  motion  is  primarily  from  one  mode,  this  can  be  simplified  by 
noting  that  [K]  -  ca*  [M]  =  0  (equation  2-68)  and  if  <|>,  is  scaled  so  that  <|)i^  [M]  <j>|  = 
1.0.  This  reduces  equation  2-69  to 


iff*  = 

dx,  ax^  *  ar,  p 


(2-76) 


after  premultiplying  both  sides  of  the  equation  by  The  frequency  is  then 
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dependent  on  both  the  variation  of  [K]  and  [M]  with  the  design  variables  whereas  the 
previous  constraints  were  solely  based  on  the  variation  of  [K]  with  respect  to  the 
design  variables. 

Representations  for  the  design  variables,  objective  function,  and  constraint 
functions  have  been  presented.  These  equations  can  now  be  used  to  calculate  the 
minimum  weight  design  for  the  vertical  tail  by  reducing  the  sizes  of  each  of  the 
structural  members  to  the  minimum  size  which  will  not  violate  any  of  the  active 
constraints. 
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3 


Vertical  Tail  Model 


3.1  Introduction 

The  mathematical  tools  developed  in  chapter  two  are  now  used  to  create  a  working 
mathematical  model  which,  theoretically,  replicates  the  structural  response  of  the  vertical 
tail  to  aerodynamic  loads.  The  mathematical  model  is  created  for  analysis  with  a  program 
called  MSC/NASTRAN  (Reference  25).  To  create  the  model,  specific  information  about 
the  design  is  necessary.  This  information  is  used  to  build  finite  element  matrices  such  as 
the  structural  stiffness  matrix  (equation  2-22)  and  the  pressure-downwash  matrix  (equation 
2-16). 

The  mathematical  model  can  be  broken  into  three  major  components.  A 
description  of  the  three  major  components  of  the  vertical  tail  model  is  given  in  the 
following  sections.  The  next  section  (3.2)  describes  the  structural  model.  The  physical 
properties,  geometry,  and  limitations  of  the  model  is  presented.  Section  3.3  explains  how 
the  optimization  problem  is  setup.  The  next  section,  3.4,  describes  the  aerodynamic  grid, 
which  generates  the  aerodynanuc  loads  that  are  used  for  the  structural  model. 

The  final  section,  3.5,  addresses  the  reasons  for  choosing  MSC/NASTRAN  as  the 
analysis  tool.  The  advantages,  as  well  as  the  disadvantages,  of  using  this  package  are 
featured  in  this  section.  Comments  on  the  performance  of  NASTRAN  in  calculating  the 
thesis  results  are  delayed  until  chapter  five. 
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3.2 


Vertical  Tail  Structure 


The  vertical  tail  is  a  highly-swept  lifting  surface.  It  is  assumed  that  the  root  of  the 
tail  is  cantilevered,  or  attached  to  a  rigid  surface.  The  tail  span  is  3  meters  and  the  total 
surface  area  is  5.4  m^  The  chord  length  tapers  from  2.4  meters  at  the  root  to  1.2  meters 
at  the  tip.  The  leading  edge  of  the  tail  is  swept  45  degrees.  The  tail  is  oriented  in  a 
vertical  fashion  with  respect  to  the  fuselage  of  the  aircraft.  Structural  support  for  the  tail 
surface  is  provided  from  a  tail  box.  The  tail  box  has  a  typical  3  spar  configuration. 


12  m 


Figure  3-1:  Votical  Tail  Dimensions. 
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The  tail  box  is  a  3  spar,  10  rib  design.  The  box  is  thickest  at  the  root  of  the 
middle  spar  (16  cm).  The  forward  and  aft  spar  are  12  cm  thick  at  the  root.  The  thickness 
of  the  wing  box  tapers  to  8  cm  at  the  tip  of  the  tail  box.  This  can  be  seen  in  figure  3-2. 
The  10  rib  sections  interconnect  the  three  spar  sections  and  are  evenly  spaced  along  the 
span  of  the  tail  box. 

The  tail  box  is  assumed  to  be  made  of  Aluminum  2017.  The  material  properties 
are  given  in  table  3-1. 


Table  3-1:  Material  Properties  of  Aluminum  2017  (T-4). 


Properties 

Values 

Modulus  of  Elasticity 

72.4(10®)  N  /  m^ 

Poisson’s  Ratio 

0.33 

Yield  Stress 

2.76(10*)  N  /  m^ 

Specific  Gravity 

2.8(10*)  kg  /  m* 

To  simplify  the  structural  analysis  only  the  tail  box  is  modeled.  The  elastic 
behavior  of  the  leading  and  trailing  edges  of  the  tail  are  not  modeled.  This  is  a 
reasonable  simplification  since  a  majority  of  the  leading  and  trailing  sections  of  the  tail 
are  directly  attached  to  the  tail  box  and  will  deform  only  slightly  with  respect  to  die  tail 
box.  The  aerodynamic  model  does  model  the  entire  tail  surface,  however,  and  the  loads 
applied  to  the  leading  and  trailing  edge  sections  are  transfened  to  the  tail  box  as  a  force- 
moment  couple.  Reducing  the  structural  model  to  just  the  tail  box  greatly  simplifies  the 
model  and  allows  the  designer  to  focus  on  the  structural  members  which  redistribute  the 
majority  of  the  applied  forces. 
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By  making  this  assumption,  however,  the  effect  of  the  aerodynamic  forces  which 
act  over  the  entire  surface  of  the  tail  must  somehow  be  transferred  to  the  tail  box.  This 
technique,  referred  to  as  splining,  was  discussed  in  section  2.3.  This  technique  transfers 
the  loads  so  that  the  moment  and  forces  produced  by  the  applied  forces  on  the  tail  box  are 
identical  to  the  forces  and  moments  produced  by  the  entire  surface  of  the  tail.  (See 
Figure  3-2) 

The  tail  box  structural  members  are  divided  into  six  groups.  The  groups  are  the 

(1)  Rib  Bars 

(2)  Spar  Bars 

(3)  Post  Bars 


Figure  3-2  :  Equivalent  transfer  of  loads  frtmi  tail  surface  to  tail  box. 
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(4)  Skin  Panels 

(5)  Spar  Panels 

(6)  Rib  Panels 

Groups  1  -  3  are  modeled  with  bar  elements  whereas  groups  4  -  6  are  modeled  with  panel 
elements.  Skin  panels  are  located  on  the  top  and  bottom  surface  of  the  tail  box  and 
provide  stiffness  against  vertical  motion  along  with  the  spar  panels.  Skin  panels  provide 
tensile  and  compressive  stiffness  whereas  the  spar  panels  provide  stiffness  with  their  shear 
strength  properties.  Resistance  to  torsional  forces  in  the  tail  box  is  provided  primarily 
from  the  shear  properties  of  the  rib  panels  which  are  positioned  perpendicular  to  the  mid¬ 
spar.  Bars  are  located  along  the  edge  of  each  of  the  panels  to  represent  the  bar  covers 
which  are  used  along  each  panel  connection. 

The  bar  elements  provide  stiffness  in  the  axial  direction  of  the  bar.  No  stiffness  is 
provided  in  any  direction  other  than  the  axis  of  the  bar.  Therefore,  forces  in  directions 
other  than  the  axis  of  the  bar  must  be  resisted  by  adjacent  elements.  The  displacement 
along  the  span  of  the  bar  element  is  assumed  to  vary  linearly.  Therefore,  from  equation 
(2-24)  and  (2-25)  the  stress  is  constant  through  each  bar  element.  The  panels  are  modeled 
with  elements  which  replicate  plane-stress  behavior.  These  are  referred  to  as  panel 
elements.  These  elements  provide  stiffness  in  the  plane  of  the  element  against  tensile, 
compressive,  and  shear  forces.  The  displacement  function  of  each  panel  element  is 
assumed  to  be  quadratic.  Therefore,  the  stresses  vary  linearly  throughout  the  element. 

The  elements  are  comiected  at  points  referred  to  as  nodal  points.  All 
displacements  are  calculated  at  these  nodal  points.  The  bar  element  requires  two  nodal 
points  per  element  Two  types  of  panel  elements  are  used.  A  common  type  is  the 


3-5 


*  =  Node  Point 


Figure  3*3  :  Description  of  location  of  each  of  the  six  groups  of  elements. 
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quadrilateral  panel  element.  This  element  requires  four  nodal  points.  Only  two  elements 
use  the  second  type  of  element  which  requires  only  three  nodal  points.  This  is  referred  to 
as  a  triangular  panel  element.  There  are  58  nodal  points  in  the  tail  box  design.  The  tail 
box  is  composed  of  80  panel  elements  and  1 14  bar  elements.  (See  rigure  3-4) 

The  stiffness  matrix  (  equation  2-22)  for  this  problem  can  be  calculated  with  the 
information  known  about  the  geometry  of  the  tail  box,  and  the  properties  of  the  bar  and 
panel  element.  The  mass  matrix  is  formed  by  assuming  the  mass  is  evenly  distributed 
throughout  the  element.  This  implies  that  the  mass  is  evenly  distributed  among  the  nodal 


—  u(x) 


Bar  EHement 


The  unknown  constants, 
Cf  are  detennined  from 
the  displacements  at  u 
and  V  J 


u(x)-Cj+C2X 


u(x.y)  -  c  1+  c  ^  +  c  73+ 
v(x.y)  -  c  5+  c  c  y^* 


Figure  3^  :  Element  Types  used  in  design  model. 
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points  of  each  element.  Since  the  dynamic  response  of  the  tail  box  is  sensitive  to  the 
mass  distribution,  additional  mass  was  added  to  various  nodal  points  to  represent  the  mass 
of  the  leading  and  trailing  edges  of  the  tail,  and  any  other  non-structural  mass  (electrical 
equipment,  paint,  etc.)  that  should  be  included  to  accurately  model  the  dynamic  behavior 
of  the  structure.  The  total  mass  added  was  34  kg. 

If  the  applied  loads  are  given,  then  enough  information  is  now  available  to  solve 
the  displacement-force  equilibrium  equation  for  the  unknown  displacements;  from  which 
the  stress,  strains,  etc.  can  be  obtained. 

3.3  Optimization  Model 

Some  additional  information  is  required  to  define  how  the  optimization  process 
will  progress.  First,  the  optimization  problem  uses  all  the  information  from  the  structural 
model.  In  addition,  the  design  variables,  objective  function,  and  constraint  functions  have 
to  be  defined.  The  objective  and  constraint  equations  must  be  declared,  but  they  do  not 
directly  affect  the  model  so  further  discussion  is  pos{poned  until  the  next  chapter.  Note, 
however,  that  the  objective  of  the  design  is  to  minimize  the  weight  and  the  only 
constraints  considered  are  structural  constraints. 

The  cross-sectional  area  of  the  bar  elements  and  the  thickness  of  the  panel 
elements  are  defined  as  the  design  variables.  This  implies  that  there  are  194  design 
variables  in  the  model.  (Total  number  of  panel  and  bar  elements)  The  complexity  and 
size  of  the  problem  increases  rapidly  as  the  number  of  design  variables  are  increased.  So, 
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to  reduce  the  problem  siz^,  design  variable  groups  are  defined  where  the  design  variables 
in  each  group  are  forced  to  be  equal  in  value  (linked  together).  Three  schemes  for  design 
variable  linking  are  used.  The  first  scheme  reduces  the  problem  to  two  design  variables. 
All  of  the  bar  elements  are  linked  together  as  one  design  variable  and  all  of  the  panel 
elements  are  linked  together  as  one  design  variable.  The  second  scheme  includes  the 
variation  of  stiffness  requirements  between  the  six  components  of  the  tail  box  defined 
earlier  by  linking  all  of  the  posts  together  as  one  design  variable,  all  of  the  spars  as  one 
design  variable,  etc.  In  all,  there  are  6  design  variables.  The  last  scheme  subdivides  each 
of  the  six  components  by  location  on  the  tail.  The  design  variables  are  given  as  follows: 


1. 

Post  -  Inboard  (Bar  Elements) , 

2. 

Post  -  MidSpan  (  ” 

"  ). 

3. 

Post  -  MidSpan  (  " 

"  ). 

4. 

Post  -  Outboard  ( " 

"  ). 

5. 

Spar  -  Inboard  ( " 

"  ), 

6. 

Spar  -  MidSpan  ( " 

"  ). 

7. 

Spar  -  MidSpan  (  " 

"  ). 

8. 

Spar  •  Outboard  (  " 

’’  ). 

9. 

Ribs  -  Inboard  (  " 

"  ), 

10. 

Ribs -MidSpan  (" 

"  ), 

11. 

Ribs  -  (Outboard  (  " 

"  ), 

12. 

Skin  Panels-Inboard 

(Panel 

Elements) 

13. 

Skin  Panels-MidSpan 

(  " 

"  ). 

14. 

Skin  Shear  Panels-Outboard  (  " 

"  ), 

15. 

Spar  Shear  Panels-Inboard  (  " 

"  ), 

16. 

Spar  Shear  Panels-MidSpan  (  " 

"  ), 

17. 

Spar  Shear  Panels-Outboard  (  " 

"  ), 

18. 

Rib  Shear  Panels-Inboard  (  " 

"  ). 

19. 

Rib  Shear  Panels-MidSpan  (  " 

"  ). 

20. 

Rib  Shear  Panels-Outboard  (  " 

"  ). 

21. 

Rib  Shear  Panels-Outboard  (  " 

"  ). 

The  location  of  the  linked  design  variables  is  given  in  figures  (3-5  through  3-10). 

Now  that  the  design  variables  are  defined,  the  optimization  code  can  use  these 
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Figure  3*5  ;  Illustration  of  the  six  and  twenty  one  design  variable  linking  scheme  for  the 
spar  (bar)  elements. 


Figure  3-6 :  Illustration  of  the  six  and  twenty  one  design  variable  linking  scheme  for  the  spar 
panel  elements. 
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6  dv  rib  group 


rib  group  3  :  21  dv 


=30 

■fclJO  ^ 


rib  group  2 


rib  group  1  :  21  dv 


Figure  3-7  ;  Illustration  of  the  six  and  twenty  (Mie  design  variable  linking  scheme  for  the  nb 
(bar)  elements. 


6  dv  rib  panel  group 


,N^/ 


rib  panel  group  3  ;  21  dv 


rib  panel  group  2  :  21  dv 


rib  panel  group  1  ;  21  dv 


Figure  3-8 :  lUustratiOT  of  the  six  and  twenty  one  design  variable  linking  scheme  for  the  rib 
panel  elements. 
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Figure  3-9 :  Illustration  of  the  six  and  twenty  one  design  variable  linking  scheme  for  the  post 
(bar)  elements. 


Figure  3-10  :  Illustration  of  the  six  and  twenty  one  design  variable  linking  scheme  for  the 
skin  panel  elements. 
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variables  in  the  objective  function  and  enforce  the  constraint  equations  on  the  design.  The 
last  item  that  needs  to  be  addressed  is  the  applied  forces.  The  optimization  code  can  only 
use  forces  which  are  independent  of  displacement.  This  becomes  a  problem  with 
aerodynamic  loads  which  are  dependent  on  structural  displacement.  To  alleviate  this 
coupling,  the  loads  used  with  the  optimization  sequence  were  the  rigid  tail  loads.  Rigid 
tail  loads  are  the  aerodynamic  loads  that  would  be  produced  if  the  wing  were  infinitely 
stiff  and  did  not  deform.  By  not  including  the  aeroelastic  influence  on  the  structural 
loading  the  accuracy  of  the  optimization  is  decreased.  The  signiricance  of  the  flexible 
aerodynamic  loads  on  a  optimal  design  based  on  rigid  tail  aerodynamics  is  addressed  in 
chapter  4. 

The  iterative  optimization  algorithm  assumes  that  an  optimal  solution  has  been 
attained  when  the  objective  function  changes  less  than  0.1  percent  between  three 
iterations.  The  constraints  are  allowed  to  be  violated  by  0.5  percent,  and  the  maximum 
number  of  allowed  iterations  before  the  program  was  terminated  was  set  at  20.  (The 
number  of  iterations  never  reached  20  during  the  analysis,  however) 

3.4  Aerodynamic  Model 

The  tail  is  assumed  to  have  a  symmetrical  surface  with  no  camber  or  twist.  This 
allows  the  tail  to  be  treated  as  a  flat  plate.  The  surface  is  positioned  at  the  mid-plane 
between  the  top  and  bottom  surface  of  the  tail  box.  To  conrigure  the  tail  in  a  manner 
which  can  use  the  doublet-lattice  method  of  section  2.2,  the  tail  is  subdivided  into  a  finite 
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Figure  3-11  :  Aerodynamic  Panel  Mesh. 

number  of  panels.  Reference  18  recommended  at  least  six  chordwise  panels  and  that  each 
of  the  panels  maintains  an  aspect  ratio  near  1.0.  This  recommendation  was  taken  at  first, 
but  it  was  found  that  by  increasing  the  numb^  of  chordwise  panels  to  seven,  and  adding 
the  appropriate  number  of  panels  in  the  span  wise  direction  to  keep  the  aspect  ratio  near 
1.0,  that  the  lift  values  varied  significantly.  Tlie  large  variation  in  lift  tapered  off  to 
negligible  differences  for  the  case  of  eight  and  then  nine  chordwise  panels,  so  the 
configuration  of  eight  chordwise  panels  was  chosen.  To  keep  the  aspect  ratio  near  one, 
ten  spanwise  panels  were  required.  This  brought  the  total  number  of  panels  to  eighty. 

The  resulting  mesh  of  aerodynamic  panels  is  shown  in  figure  3-11.  Each  of  the  panels 
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has  a  doublet  line  at  the  quarter  chord  line  and  a  control  point  at  the  mid-span  of  the 
three-quarter  chord  line,  resulting  in  80  control  points  from  which  to  solve  the  downwash- 
pressure  equation  (equation  (2-17). 

Each  of  these  panels  is  treated  as  a  rigid  flat  plate  with  freedom  to  displace 
vertically  and  rotate  about  the  control  point.  If  an  initial  angle  of  attack  is  given,  then  the 
entire  surface  is  rotated  about  a  given  point  at  the  root  of  the  tail  and  the  iititial  downwash 
for  each  panel  would  be  the  slope  of  the  individual  panel  with  respect  to  the  x-coordinate 
direction  (see  equation  2-4).  As  the  tail  box  deforms  under  the  applied  loads,  these 
stmctural  deformations  are  then  converted  to  appropriate  aerodynamic  deformations  at  the 
aerodynamic  control  points  through  the  splining  matrix,  and  appropriate  changes  are  made 
to  the  relative  vertical  displacement  and  rotation  of  each  panel. 

3.5  Comments  on  MSC/NASTRAN 

The  theory  discussed  in  chapter  3  is  incorporated  into  a  program  marketed  by  the 
MacNeal-Schwendler  Corporation  (Reference  18).  This  program  is  called 
MSC/NASTRAN  (Version  66-A).  The  program  was  used  to  produce  the  results  given  in 
this  thesis.  This  program  was  chosen  because  of  it’s  structural  analysis,  aeroelastic 
analysis,  and  optimal  design  capabilities.  To  accomplish  each  of  these  tasks,  at  least  one 
of  three  major  data  decks  had  to  be  used.  The  first  deck  contains  all  of  the  structural 
information.  For  example,  structural  grid  point  information,  element  description,  element 
connectivity,  material  type,  etc.  The  second  deck  contained  all  of  the  structural 
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Figure  3*12  :  The  Three  NASTRAN  Modules  utilized  for  the  analysis  of  the  vertical  tail. 


information  and  additionally  aerodynamic  information  such  as  flight  conditions, 
aerodynamic  mesh  description,  planes  of  symmetry,  etc.  The  last  deck  described  the 
optimization  design  model.  Therefore,  the  objective  function,  constraints,  and  design 
variables  had  to  be  defined.  The  optimization  is  based  on  the  structural  model  so  the 
structural  information  had  to  be  included  in  the  model.  The  modular  structure  of 
NASTRAN  made  the  individual  analysis  siQq)ler,  which  allowed  the  author  to  track  the 
flow  of  information  more  effectively.  Examples  for  each  of  the  data  decks  is  provided  in 
appendix  B,  C,  and  D. 

The  loads  used  in  the  aeroelastic  analysis  are  based  on  doublet-lattice 
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aerodynamics  (1:279).  Therefore,  the  flutter  analysis  uses  a  paneling  method  based  on 
unsteady,  potential  aerodynamic  flow.  The  loads  calculated  reflect  the  change  in 
aerodynamic  properties  due  to  the  flexible  nature  of  the  structure.  This  applies  to  both  the 
static  and  dynamic  response. 

The  downside  of  NASTRAN  was  that  their  is  no  "bridge"  between  optimal  design 
and  aerodynamic  analysis.  For  example,  the  optimal  design  is  based  on  the  aerodynamic 
loads  of  a  rigid  tail.  These  loads  are  extracted  from  the  output  of  the  static  aeroelastic 
analysis  by  segregating  the  loads  due  to  flexibility  from  the  loads  due  to  initial 
conditions.and  then  inserted  into  the  optimization  analysis.  This  deficiency  stems  from 
the  inability  of  the  NASTRAN  code  to  update  the  change  in  the  aeroelastic  loads 
structural  deformations  during  the  optimization  procedure.  Another  restriction  or  i 

NASTRAN  is  that  the  constraint  functions  are  restricted  to  structural  properties.  Thex 
no  opportunity  to  constrain  the  values  of  the  stability  derivatives,  local  angle  of  attack, 
etc. 

A  program  without  these  deficiencies  would  enhance  the  accuracy  of  the  analysis, 
but  unfortunately  a  program  with  these  capabilities  was  not  available  to  analyze  a 
cantilevered  lifting  surface. 
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4 


Results 


4.1  Introduction 

This  chapter  presents  the  results  from  the  analysis  of  the  vertical  tail.  The  analysis 
focussed  on  the  influence  that  design  parameters  have  on  the  optimal  size  and  weight  of  the 
tail  box.  These  parameters  include  the  choices  of  design  variable  linking  scheme,  critical 
design  load,  and  design  constraints.  The  influence  of  each  factor  is  presented  in  the  first 
three  sections. 

The  aeroelastic  properties  of  each  of  the  designs  are  then  evaluated.  Comparisons  of 
the  aerodynamic  properties  of  each  of  the  designs  are  discussed.  Section  4.5  focuses  on  the 
results  of  the  static  aeroelastic  analysis  while  section  4.6  presents  the  flutter  analysis  results. 

A  discussion  and  comparison  of  the  optimization  and  aerodynamic  results  will  follow  in 
chapter  five. 

4.2  Effect  of  Design  Variable  Linking  on  Objective  Function 

Three  design  linking  schemes  are  considoed  to  illustrate  the  influence  of  design 
variable  linking  on  the  results  of  an  optimization  analysis.  Design  variable  linking  is  valuable 
because  of  the  savings  in  computer  time.  The  three  linking  schemes  described  in  section  3.3 
are  used.  Recall  the  first  scheme  subdivided  the  model  into  2  groups,  the  bar  elements  make 
up  one  design  variable  group  and  the  panel  elements  matfe  up  the  other  group.  The  second 
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scheme  further  subdivided  the  bar  elements  and  panel  elements  into  3  design  variable  groups 
respectively,  (six  groups  total)  The  final  scheme  was  yet  a  further  subdivision  of  the  second 
scheme,  and  is  composed  of  21  design  variable  groups.  These  groups  will  be  periodically 
referred  to  as  design  variables  themselves.  (See  Figures  3.3  -  3.8) 

Each  linking  scheme  is  compared  by  observing  the  resulting  objective  weight,  design 
variable  size,  and  element  stress  values.  Only  one  loading  condition  is  considered  for  each 
case  and  the  only  constraint  enforced  on  the  design  is  maximum  stress  allowed.  The  tail  load 
is  calculated  at  Mach  number  0.4  (sea  level)  and  30  degrees  angle  of  attack  (steady  level 
flight). 

The  optimal  weight  for  each  case  was  obtained  from  the  optimization  output.  The 
heaviest  design  was  the  two  design  variable  case  which  had  an  optimal  weight  of  73.5  kg. 

By  increasing  the  number  of  design  variables  to  six,  the  weight  of  the  design  dropped 
considerably  to  44.3  kg.  The  21  design  variable  case  subdivided  each  of  the  six  design 
variables  into  three  or  four  subgroups.  By  using  the  21  design  variable  scheme,  the  optimal 
weight  calculated  dropped  even  further  to  38.2  kg. 

To  illustrate  where  the  structural  mass  is  being  used  to  create  a  feasible  design,  the 
values  for  each  of  the  design  variables  (member  sizes)  are  compared  graphically.  Since  the  6 
and  21  design  variable  case  are  subdivisions  of  the  2  design  variable  case,  comparisons  of 
member  sizes  are  easily  shown.  These  comparisons  are  shown  in  frgures  4-1  through  4-6. 

The  stress  values  along  three  selected  tail  box  regions  are  also  graphically  depicted  to 
illustrate  the  dependence  of  element  stress  values  on  the  design  variable  scheme  chosen. 

Three  groups  of  elements  were  chosen.  The  first  group  of  elements  is  the  middle  spar 
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of  Panel  (10(-3)m) 


Figure  4-1 :  Optimal  design  variable  size  •  cross-sectional  area  of  spar  elements. 


Figure  4-2  :  Optimal  design  variable  size  -  thickness  of  spar  panel  elements. 


4-4  :  Optimal  design  variable  size  •  thickness  of  rib  panel  elements. 
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Figure  4-5  :  Optimal  design  variable  size  -  cross-sectional  area  of  post  elements. 
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Figure  4-6  :  Optimal  design  variable  size  -  thickness  of  skin  panel  elements. 
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1-7  :  Von  Mises  stress  values  -  mid  spar  panel  elements. 
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Figure  4-8  ;  Von  Mises  stress  values  -  forward  rib  panel  elements. 
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Element  Number 

Figure  4-9  :  Von  Mises  stress  values  -  forward  skin  panel  elements. 

section.  The  Von  Mises  stress  for  the  string  of  spar  panel  elements  from  the  root  chord  to 
the  tip  chord  (left  to  right)  are  depicted  in  Figure  4-7.  Note  that  the  stress  values  given  are 
the  stresses  at  the  geometric  centroid  of  the  element.  The  second  group  is  the  forward  section 
of  rib  panel  elements.  The  elements  are  located  along  the  span  of  the  tail  box  from  root  to 
tip  also  (Figure  4-8).  The  last  group  is  the  forward  section  of  the  skin  panels,  these  elements 
also  are  presented  in  Figure  4-9  ordered  from  root  to  tip  section  of  the  tail  box. 

4.3  Effect  of  Design  Load  Condition  on  Minimum  Weight 

A  relationship  between  the  design  load  and  the  required  structural  weight  to  sustain 
the  load  could  be  valuable  in  determining  the  weight  penalty  for  an  increased  tail  design  loaH 
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This  relationship  is  examined  by  optimizing  the  tail  for  a  variety  of  possible  design  load 
conditions. 

Again,  the  tail  loads  are  based  on  rigid  tail  aerodynamics  and  the  optimization  is  based 
on  the  linking  scheme  of  six  design  variables.  The  six  design  variable  case  is  chosen  since  it 
represents  a  considerable  improvement  of  the  objective  function  value  over  the  two  design 
variable  case  yet  does  not  require  near  the  computer  time  that  the  21  design  variable  case 
does.  The  structure  is  optimized  for  minimum  weight  and  the  design  variables  are 
constrained  by  maximum  stress  limits  and  minimum  gauge  limits.  The  minimum  gauge  limits 
are  enforced  to  examine  how  the  minimum  sizes  effect  the  optimal  weight  values. 

The  load  cases  were  determined  by  providing  the  Mach  number  and  angle  of  attack 
parameters.  All  cases  assume  the  tail  is  at  sea  level  in  steady  level  flight.  The  cases  vary 
from  a  Mach  number  of  0.2  to  a  Mach  number  of  0.8.  The  peak  design  angle  of  attack  was 
40  degrees  for  M  =  0.2,  but  tapered  to  6  degrees  for  M  =  0.8.  It  is  recognized  that  the  results 
at  high  angles  of  attack  are  tainted  since  the  aerodynamic  assumption  that  the  surface  is  at  a 
small  angle  of  attack  is  violated.  This  error  is  accepted  and  noted  for  this  analysis.  The 
design  weight  varied  from  16.2  kg  (all  members  minimum  gauge)  for  the  M=0.2  a  =  10  deg, 
to  58.9  kg  for  the  case  of  M=0.6,  a=15  deg.  The  variation  in  weight  is  illustrated  in  figures 
4-10  and  4-11.  Figure  4-10  shows  the  linear  variation  of  design  weight  with  angle  of  attack 
for  a  given  mach  number.  Figure  4-11  better  depicts  the  linear  relationship  between  tail 
aerodynamic  load  and  design  weight  The  lines  are  not  co-linear  but  definitely  show  a 
general  linear  relationship  between  design  load  and  optimal  design  weight.  The  cause  for  the 
large  deviation  along  the  M=0.4  line  cannot  be  explained. 
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Figure  4>10  :  Optimal  structural  weight  versus  design  angle  of  attack  for  various  flight 
conditions  at  sea  level. 
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Figure  4-11 :  Optimal  weight  of  design  compared  to  the  applied  aerodynamic  load  used  for 
each  of  the  optimizatitm  cases. 
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4.4  Effect  of  Constraints  on  Final  Weight 


To  demonstrate  the  influence  that  design  constraints  have  on  an  optimal  design,  the 
tail  is  optimized  for  three  different  cases.  The  first  case  constrained  the  values  for  the 
element  sizes  so  that  the  maximum  stress  does  not  exceed  the  material  yield  stress.  The 
second  case  adds  an  additional  restriction  that  the  member  sizes  can  not  decrease  below  a 
minimum  value.  The  last  case  involves  an  attempt  to  drop  the  tip  acceleration  on  the  tail  by 
constraining  the  tip  deflection.  Only  one  loading  condition  is  used  for  each  case;  M  =  0.4  at 
an  angle  of  attack  of  30  degrees  (sea  level,  steady  level  flight,  and  rigid  wing  aerodynamics). 
In  addition,  to  illuminate  any  spanwise  differences  between  each  of  the  designs,  the  linking 
scheme  of  21  design  variables  is  used. 

The  first  case  enforces  internal  stress  limitations  of  the  material.  For  the  bar  element, 
each  bar  cross-section  is  not  allowed  to  exceed  the  yield  stress  of  2.76  (10*)  nW.  The  panel 
element  thickness  values  are  not  allowed  to  drop  below  a  value  which  would  cause  the  Von 
Mises’  stress  to  exceed  2.76  (10*)  N/m^.  The  stress  constraint  (sc)  restricted  the  optimal 
weight  from  dropping  below  38.2  kg.  After  completing  the  optimization,  the  stress  values 
indeed  remained  below  the  maximum  allowed.  (See  Hgure  4-18) 

The  second  constraint  case  enforced,  in  addition  to  maximum  stress,  a  minimum  gauge 
(mg)  constraint  to  model  a  manufacturing  limitation  which  might  limit  the  size  of  the  panel  or 
bar  firom  dropping  below  a  given  value.  The  minimum  gauge  (mg)  values  chosen  for  this 
analysis  are  0.1  cm^  (1.  10*^  m^)  cross-section  for  the  bar  elements  and  1.0  mm  (1.  10'^  m) 
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Figure  4-15  :  Optimal  design  variable  size  •  thickness  of  rib  panel  elemrats. 
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Figure  4-16  :  Optimal  design  variable  size  -  cross-sectional  area  of  post  elements. 
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Figure  4-17  :  Optimal  design  variable  size  -  thickness  of  skin  panel  elements. 
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thickness  for  the  panel  elements.  This  constraint  increased  the  size  of  the  post  elements 
(Figure  4-13)  and  skin  panels  (Figure  4-16).  Because  of  the  increase  in  size,  the  stress  level 
in  some  of  the  elements  which  are  near  the  minimum  limit  in  the  stress  constrained  case 
decreased,  as  is  seen  for  the  mid-spar  elements  of  figure  4-18.  The  weight  penalty  for  this 
constraint  is  minor,  however,  as  the  weight  only  increased  to  38.4  kg. 

The  third  constraint  case  additionally  restricted  the  tail  tip  (nodal  point  28)  from 
deflecting  more  than  10  cm.  For  the  minimum  gauge  case  (case  two),  the  deflection  was  19.6 
cm.  Why  would  a  reduction  in  tip  displacement  be  necessary?  Electronic  equipment  which 
is  often  placed  at  the  tail  tip  is  sensitive  to  accelerations.  If  a  piece  of  equipment 
(navigational  or  electronic  countermeasure  for  example)  must  not  ’see’  accelerations  greater 


Figure  4*18  :  Comparison  of  stress  values  in  identical  skin  panel  dements  for  stress 
constrained  case  (sc)  and  stress  constrained/  minimum  gauge  (mg)  case. 
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than  a  given  value,  then  the  structure  must  be  redesigned  to  insure  that  the  maximum  allowed 
acceleration  is  not  exceeded.  The  accelerations  at  the  tail  tip  are  assumed  to  be  proportional 
to  the  displacement  at  the  tip.  (  acceleration  =  (firequency)^  x  vertical  displacement )  This  is 
true  for  harmonic  motion.  If  the  tail  can  be  redesigned  to  reduce  this  deflection  then  the 
maximum  acceleration  should  decrease  linearly.  This  assumption  turned  out  to  be  incorrect 
because  as  the  wing  is  stiffened  to  reduce  the  tip  deflection,  the  added  stiffness  also  increases 
the  first  bending  frequency.  The  displacement  constraint  dropped  the  deflection  from  19.6  cm 
to  10.0  cm,  however,  this  constraint  increased  the  tint  frequency  (a  first  bending  mode)  from 
11.30  Hz  (sc/mg)  to  16.91  Hz.  The  result  was  a  net  gain  in  tip  acceleration  from  988.03  m/s^ 
for  the  minimum  gauge  case  to  1 128.9  m/s^  for  the  displacement  constraint  (dc)  case. 
Therefore,  this  type  of  constraint  was  not  successful  in  reducing  the  tip  acceleration  and, 
therefore,  a  relationship  between  the  tip  acceleration  and  the  design  variables  will  need  to  be 
developed  to  reduce  the  tip  acceleration. 

It  is  interesting  to  note  where  the  optimization  code  increased  the  sizes  of  the 
members  to  reduce  the  tip  deflection.  The  members  which  were  affected  most  were  the  spar 
panels,  rib  panels,  and  skin  panels.  The  displacement  constraint  (dc)  also  increased  the 
optimal  weight  to  40.5  kg. 

Now  after  examining  the  structural  behavior  of  these  different  designs,  the  data  which 
is  most  important  to  the  performance  of  the  aircraft  is  examined  next.  This  data  is  the 
aeroelastic  response,  static  and  dynamic,  of  the  tail  to  given  aerodynamic  loads. 
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4.5  The  Static  Aeroelastic  Properties  of  Constrained  Designs. 


The  aerodynamic  properties  of  a  lifting  surface  are  closely  related  to  the  stmctural 
response  to  applied  aerodynamic  loads.  To  illustrate  this,  the  structiual  displacements  of  each 
of  the  constrained  designs  of  section  4.4  are  related  to  three  basic  static  stability  derivatives 
’  '^^®re  a  is  the  angle  of  attack  of  the  tail.  The  symbol  is  the  lift 

coefficient  which  relates  angle  of  attack  to  the  lift  on  the  tail  (  Lift  =  q  S  C,  a).  The 

symbols,  and  C„  ,  relate  the  pitch  moment  about  the  V4  chord  line  and  roll  moment 

about  the  tail  root  span  to  angle  of  attack  (  Pitch  Moment  =  q  S  c  C*.  a  and  Roll  Moment 

=  q  S  b  Cff^  a),  where  c  is  the  root  chord  length  and  b  is  the  full  span  of  the  tail.  Recall 

that  the  structure  was  optimized  with  respect  to  the  rigid  tail  loads.  The  true  aerodynamic 
loads  (tail  flexibility  accounted  for)  differ  from  the  loads  used  to  optimize  the  structural 
members.  The  difference  in  loads  affect  the  displacement  and  member  stresses  of  each 
design. 


TABLE  4-1  :  Aerodynamic  Stability  Derivatives. 


1  Case 

Load(N) 

1  Rigid  Tail 

3.269 

-1.002 

0.7386 

104892.1 

1  Stress  Constraint 

3.498 

-1.1168 

0.7992 

112229.5 

Minimum  Gauge 

3.269 

-1.005 

0.7408 

104874.1 

Displacement  Constraint 

3.370 

-1.0514 

0.7648 

108129.1 
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An  aeroelastic  analysis  is  conducted  for  the  three  constraint  design  cases  of  the 
previous  section,  (stress  constrained  (sc),  minimum  gauge  (mg),  and  displacement  constraint 
(dc)  designs).  The  flexible  static  stability  derivatives  for  each  of  these  cases  are  presented  in 
table  4-1.  For  each  of  the  cases,  however,  the  change  was  less  than  7  percent.  The  stress 
constrained  case  exhibited  the  greatest  variation  of  lift  coefficient  compared  to  the  rigid  tail 
aerodynamics  lift  coefficient  (increase  from  3.269  to  3.498)  The  moment  and  rolling 
coefficients  increased  proportionaly  more  for  each  case  as  compared  to  the  overall  lift, 
suggesting  that  the  majority  of  the  added  lift  came  from  the  outboard  section  of  the  tail. 

Since  the  aerodynamic  loads  change,  it  is  expected  that  the  structural  properties  also 
change.  The  displacement  of  the  tail  due  to  the  flexible  aerodynamic  loads  changed  very 
little  as  compared  to  the  rigid  load  displacements  for  each  of  the  cases.  As  an  example, 
figures  4-19  through  4-24  show  the  small  difference  in  displacement  for  the  different  design 
cases.  For  the  stress  constrained  case,  since  the  lift  coefficient  increased  when  flexibility  was 
accounted  for,  it  is  intuitive  that  the  total  lift  increased  and,  therefore,  the  stresses  must  have 
increased  in  the  tail  box.  The  result  of  this  increase  is  illustrated  by  comparing  the  rigid  load 
stresses  against  the  flexible  load  stresses  for  the  rib  panel  elements  and  the  skin  panel 
elements  in  figures  4-25  and  4-26. 

The  static  aeroelastic  properties  of  the  tail  are  important,  but  equally  important  is  the 
maximum  speed  the  tail  can  operate  without  encountering  dynamic  instabilities  or  flutter. 

The  next  section  addresses  the  flutter  characteristics  of  each  of  the  constraint  designs. 
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Figure  4-19  :  The  stress  constrained  (sc)  case  displacements  at  nodal  points  along  the  three 
spars  due  to  rigid  tail  loads. 


Figure  4-20 :  the  stress  constrained  case  displacements  at  nodal  points  along  the  three  spars 
due  to  aeroelastically  corrected  tail  loads. 
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Figure  4-21 :  The  minimum  gauge  case  displacements  at  nodal  points  along  the  three  spars 
due  to  rigid  tail  loads. 


Figure  4-22  :  The  minimum  gauge  case  displacements  at  nodal  points  along  the  three  spars 
due  to  aeroelastically  corrected  tail  loads. 
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Figure  4*23  :  The  displacement  constraint  case  -  displacements  at  nodal  points  along  the 
three  spars  due  to  rigid  tail  loads. 


Figure  4-24  :  The  displacement  constraint  case  -  displacements  at  nodal  points  along  the 
three  spars  due  to  aeroelastically  corrected  tail  loads. 
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Figure  4*25:  Comparison  of  Von  Mises  stress  values  in  mid  spar  panels  with  rigid 
aerodynamic  loads  and  aeroelastically  corrected  flexible  Ir 
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4.5  Effect  of  Design  Constraints  on  Flutter  Speed 


The  dynamic  characteristics  of  each  of  the  three  design  are  now  examined.  Since 
dynamic  response  is  closely  tied  to  the  structural  mass,  additional  mass  is  evenly  distributed 
to  selected  tail  box  nodal  points  to  account  for  the  sections  of  the  tail  which  are  not  included 
in  the  finite  element  model  and  to  include  the  mass  of  any  equipment  which  may  be  stored  in 
the  tail.  The  location  of  the  added  mass  is  shown  in  figure  4-27.  The  total  mass  added  was 
34  kg.  (  2  kg  at  each  node)  To  further  examine  the  dependence  of  the  dynamic  response  to 
the  mass  distribution,  the  location  of  the  concentrated  mass  was  varied  in  the  minimum  gauge 
design.  The  modified  mass  distributions  are  noted  as  (1)  the  mass  forward  and  (2)  the  mass 
aft  case.  (See  Figures  4-28  and  4-29) 

The  first  four  natural  frequencies  and  mode  shapes  were  calculated  for  each  design. 
The  frequencies  are  given  in  table  4-2  and  the  mode  shapes  are  given  in  appendix  A.  The 
first  two  modes  of  each  design  case  are  predominantly  bending  modes  for  each  of  the  cases 
except  for  the  stress  constrained  case  whose  second  mode  is  predominantly  a  torsion  mode. 
The  dominant  bending  modes  are  expected  for  a  swept  tail  (wing)  (Reference  8)  The  first 
frequency  varies  from  10.2  Hz  for  the  stress  constrained  design  tt)  16.9  Hz  for  the 
displacement  constraint  design.  This  compares  well  with  the  frequency  of  the  F- 16  vertical 
tail  (  »  12.7  Hz)  and  the  F-15  fighter  aircraft  (  »  9.1  Hz).  The  higher  frequencies  of  the 
vertical  tail  model  grow  much  more  rapidly  than  the  values  for  the  two  fighter  aircraft  used  as 
a  benchmaik,  however.  For  example,  the  second  frequency  of  the  stress  constrained  design 
(sc)  is  41.28  Hz  whereas  the  F-16  fighter  aircraft  vertical  tail  has  a  second  frequency  of  »  32 
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Figure  4>27:  Normal  distribution  of  concentrated  mass. 


Figure  4-28:  Mass  forward  distribution  of  concentrated  mass. 
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Figure  4-29:  Mass  aft  distribution  of  concentrated  mass. 


Hz. 

Table  4-2  :  Natural  Frequencies  of  the  Different  Designs. 


Design  Condition 

1st  Frequency 
(Hz) 

2nd 

Frequency 

(Hz) 

3rd 

Frequency 

(Hz) 

4th  Frequency 
(Hz) 

Stress  Constrained 
(sc) 

10.20 

41.28 

53.94 

89.49 

Minimum  Gauge 

(mg) 

11.30 

43.01 

55.23 

97.70 

Node  28  <  10  cm 
(dc) 

16.91 

50.55 

56.79 

106.26 

Weight  Forward  (mg) 

12.65 

42.78 

75.41 

90.58 

Weight  Aft  (mg) 

13.59 

43.42 

72.80 

100.8 

The  fundamental  flutter  speed  in  each  of  the  design  cases  is  predominantly  affected  by 
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the  second  mode  (bending  mode).  Equation  2-47  was  used  to  determine  the  flutter  speed. 

The  aerodynamic  matrix,  [A],  is  calculated  for  a  series  of  reduced  frequencies  at  Mach 
number  of  0.8.  (sea  level)  The  reduced  frequency,  k,  was  varied  to  determine  the  value  at 
which  the  artificial  damping  term,  g,  crossed  zero  along  with  the  corresponding  velocity. 
Graphs  of  g  for  different  values  of  velocity  and  the  forced  frequencies  of  vibration  for 
different  values  of  k  are  shown  in  figures  4-30  through  4-39.  The  lowest  value  of  velocity 
for  which  g  was  zero  was  defined  to  be  the  flutter  speed.  A  linear  interpolation  was  used 
between  data  points  to  determine  the  flutter  speed.  The  smallest  flutter  speed  calculated  was 
for  the  stress  constraint  design  (417.17  m/s).  The  minimum  gauge,  mass  forward  case 
recorded  the  highest  flutter  speed  of  737,14  m/s.  (See  Table  4-3)  These  flutter  speeds  are 
well  above  the  speed  of  sound  (  M  =  1.0  ).  To  solve  equation  2-47  exactly,  a  Mach  number  - 
velocity  match  should  be  found.  In  other  words,  the  Mach  number  at  which  the  aerodynamic 
matrices  were  calculated  should  match  the  Mach  numbo*  at  which  flutter  instability  is 
predicted.  It  is  assumed  in  this  analysis  that  the  accuracy  of  the  doublet-lattice  method  falls 
off  sharply  as  the  Mach  numbers  become  transonic  and  the  calculations  for  flutter  speed  loose 
credibility  above  M  =  0.8.  Therefore,  no  attempt  was  made  to  match  the  flutter  speed  to  the 
aerodynamic  Mach  number.  Instead,  it  is  assumed  that  flutter  is  not  important  to  the  defined 
mission  of  the  vertical  tail.  To  examine  the  flutter  mechanism  for  each  design,  velociQ^- 
damping  (v-g),  reduced  frequency  -  frequency  (k-co)  plots  are  given  for  each  design  in  figures 
4-30  to  4-39. 
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Table  4-3  :  Flutter  Speeds  for  the  Different  Designs. 

“  ,  — .  .  1 

Case  Description 

Flutter  Speed 
(m/s) 

Flutter 

Frequency 

(Hz) 

Stress  Constraint  (sc) 

417.17 

22.5 

Minimum  Gauge  (mg) 

500.17 

18.73 

Node  28  <  10  cm  (dc) 

467.64 

22.85 

Mass  Forward  (mg) 

737.14 

45.40 

Mass  Aft  (mg) 

479.99 

25.81 

Through  this  chapter  the  characteristics  and  properties  of  several  designs  have  been 
presented.  A  minimum  weight  structure  was  obtained  for  each  design  case  by  reducing  the 
sizes  of  the  structural  members.  Constraints  on  the  design  limited  the  minimum  member  size 
which  would  be  allowed.  The  structural  analysis  was  based  on  a  computer  model  which 
assumed  that  the  tail  structure  was  cantilevered  at  the  base  and  deformed  elastically  in 
response  to  the  applied  aerodynamic  loads.  The  aerodynamic  loads  were  calculated  based  on 
the  doublet-lattice  harmonic  flow  theory. 

Several  designs  have  been  evaluated  and  subtle  differences  have  been  found  between 
each  of  the  designs.  The  next  chapter  will  now  compare  and  contrast  these  differences  and 
identify  the  properties  which  would  be  used  when  considering  which  design  would  integrate 
best  into  an  aircraft  systeoL 
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Figure  4-30  :  Velocity  -  artificial  damping  (V-g)  curve  for  the  first  four  modes  of  the  stress 
constrained  (sc)  case. 
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Figure  4-31:  Reduced  frequency  •  frequency  (K-tO)  curves  for  the  first  4  modes  in  the  stress 
constrained  (sc)  case. 


4-27 


0  200  400  BOO  000  1000 

Velocity  (ni/e> 

Figure  4-32:  Velocity  •  artiflcial  damping  (V-g)  curves  for  the  first  4  modes  in  the  minimum 
gauge  (mg)  case. 
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Figure  4-33:  Reduced  frequency  •  frequency  (k-TD)  curves  for  the  first  4  modes  in  the 
minimum  gauge  constraint  (mg)  case. 


t-28 


0  200  400  BOO  BOO  1000 

VriodQr 

Figure  4>34:  Velocity  -  artiflcial  damping  (V-g)  curves  for  the  first  4  modes  in  the 
displacement  constraint  case. 
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Figure  4-35:  Reduced  frequency  *  frequency  (k-CJ)  curves  for  the  first  4  modes  in  the 
displacement  omstraint  (dc)  case. 
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Figure  4-36:  Velocity  •  artiflcal  damping  (V-g)  curves  for  the  first  4  modes  in  the  mass 
forward  (mg)  case. 
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Figure  4-37:  Reduced  frequency  -  frequency  (k-t3)  curves  for  the  first  4  modes  in  the  mass 
frarward  (mg)  case. 
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Figure  4-38:  Velocity  -  artificial  damping  (V-g)  curves  for  the  first  4  modes  in  the  mass  aft 
(mg)  case. 
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Figure  4-39:  Reduced  frequency  •  frequency  (k-GJ)  curves  for  the  first  4  modes  in  the  mass 
aft  (mg)  case. 
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5.  Discussion  and  Conclusions 


5.1  Introduction 

The  previous  three  chapters  have  been  devoted  to  developing  and  applying 
structural  optimization  to  the  design  of  the  vertical  tail.  A  considerable  amount  of 
theoretical  work  had  to  be  covered  in  chapter  two  to  describe  the  process  of  numeric 
optimization.  Chapter  three  showed  how  the  paper  design  was  transformed  into  the 
finite  element  analysis  model  and  then  modified  to  incorporate  the  optimization 
techniques  of  chapter  two.  After  considering  the  manhours  required  to  develop  these 
models,  and  the  large  amounts  of  computer  time  required  for  the  analysis  and 
reanalysis,  the  question  has  to  be  raised  whether  optimization  is  worth  the  effort  The 
results  of  chapter  four  clearly  demonstrate  that  structural  optimization  is  worth  the 
initial  cost.  This  belief  is  echoed  in  many  papers  (39:56  and  34:16). 

Altiiough  this  thesis  is  limited  to  structural  optimization,  the  development  of 
the  aeroelastic  theory  reveals  the  opportunity  to  expand  the  optimization  techniques  so 
as  to  include  aeroelastic  tailoring.  This  capability  is  in  development  today  with 
programs  such  as  ASTROS  (30:1021).  More  discussion  on  further  research  is  saved 
for  the  final  section  of  this  chapter.  The  attention  is  first  directed  toward  the 
significant  results  of  chapter  4. 
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5.2  Structural  Optimization 


Recall  that  optimization  works  by  calculating  a  value  or  series  of  values  for  the 
design  variables  which  minimize  the  objective  of  the  design.  The  objective  could  be 
anything  that  is  a  function  of  the  design  variables.  The  objective  chosen  for  the  tail 
design  was  to  minimize  the  structural  weight.  Because  of  the  matrix  formulation 
chosen  to  automate  the  optimization  process  (FEM),  the  size  of  the  matrices  grew 
considerably  as  the  number  of  design  variables  increased.  Therefore,  the  computer 
time  required  to  solve  the  problem  also  increased  rapidly  with  matrix  size.  To  keep 
the  problem  size  manageable,  the  idea  of  design  variable  linking  was  introduced  in 
chapter  3.  Three  different  linking  schemes  were  used  on  the  tail  design  to  document 
the  effect  that  linking  had  on  the  optimal  weight.  By  linking  the  elements  together, 
the  number  of  design  variables  could  be  reduced  from  197  design  variables  to  (1)2 
design  variables,  (2)  6  design  variables,  and  (3)  21  design  variables. 

There  was  a  large  difference  in  optimal  weight  between  the  two  design  variable 
and  six  design  variable  linking  scheme  (73  kg  ->  44  kg).  The  large  difference  in 
weight  was  due  to,  for  the  bar  elements,  the  highly  stressed  spar  elements  at  the  root 
of  the  tail  which  prevented  the  optimization  algorithm  from  reducing  the  size  of  the 
lesser  stressed  bars.  When  the  spar  elements  were  subdivided  into  their  own  group, 
the  other  elements  then  were  allowed  to  reduce  in  size.  (See  Figures  4-1  through  4-3) 
The  obvious  question  to  ask  after  examining  these  figures  is,  why  is  not  one  of  the 
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sizes  in  the  six  design  variable  scheme  equal  to  the  size  of  the  two  design  variable 
scheme.  The  answer  lies  in  the  complex  stress  distribution  throughout  the  tail  box.  If 
any  of  the  element  sizes  in  the  design  is  changed  then  the  stress  distribution  can 
change  drastically. 

A  more  subtle  decrease  in  the  optimal  weight  occurred  between  the  six  and 
twenty-one  design  variable  linking  scheme  (44  kg  ->  38  kg)  suggesting  that  the 
required  sizes  differ  little  at  the  wing  tip  versus  the  wing  root  The  small  change 
between  the  six  and  twenty-one  design  variable  case  was  surprising  since  it  was 
believed  that  the  spanwise  variation  in  internal  forces  would  allow  a  large  drop  in 
required  structural  weight.  It  might  then  be  concluded  that  if  minor  variations  were 
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Figure  5>1  :  Qualitative  view  of  the  effect  of  linking  together  design  variables. 
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made  to  the  design,  the  six  design  variable  linking  scheme  would  give  accurate  results 
for  a  fraction  of  computer  time. 


The  effectiveness  of  the  optimization  was  verified  by  checking  the  stress  values 
of  different  sections  of  the  tail  box.  The  direct  relationship  between  element  stress 
values  and  element  weight  allows  this  to  be  true.  The  fully  optimized  tail  stmcture 
with  only  stress  constraints  would  have  stress  values  at  the  3deld  stress  values.  This 
could  not  be  realized  because  of  the  design  variable  linking  that  was  utilized. 


Table  5-1:  Stress  value  of  the  element  with  maximum  stress  from  each  of  the  twenty  one 
design  variable  groups.  The  design  was  stress  constrained  only. 


Element 

Number 
(Bar  Element) 

Axial  Stress 

(10’  N/m*) 

Element 

Number 

(Panel 

Element) 

Von  Mises 
Stress 

il(f  N/m*) 

Spar  4 

1.42 

Skin  47 

1.94 

Spar  9 

1.15 

Skin  54 

2.74 

Spar  IS 

0.64 

Skin  57 

2.22 

Spar  21 

0.86 

Spar  232 

2.26 

Rib  27 

0.82 

Spar  235 

1.96 

Rib  35 

0.20 

Spar  241 

2.04 

Rib  39 

0.23 

Spar  247 

2.03 

Post  201 

0.10 

Rib  254 

0.75 

Post  206 

0.04 

Rib  263 

0.15 

Post  215 

0.05 

Rib  265 

0.15 

Post  226 

0.06 

However,  at  least  one  element  out  of  every  design  variable  group  should  be  fully 
stressed.  For  reasons  unknown  the  algorithm  would  halt  and  declare  that  the  design 
had  been  optimized  before  one  of  the  elements  inside  of  each  group  would  reach 
maximum  stress  of  2.76(10*)  N/m^  (see  figure  5-1).  Because  of  the  large  range  of 
numbers  used  during  the  analysis,  it  is  suspected  that  loss  in  significant  figures  during 
calculations  may  be  one  source  of  error.  This  belief  is  supported  by  the  fact  that  the 
bar  elements  did  not  optimize  as  well  as  the  panel  elements.  The  bar  elements  had 
values  three  orders  of  magnitude  less  than  the  panel  elements. 

A  relationship  was  developed  between  the  optimal  weight  and  the  design  load 
by  optimizing  the  structure  for  a  variety  of  design  loads.  The  linear  relationship  is 
shown  in  tigure  4-11.  The  linear  relationship  was  expected  since  the  stiffness 
properties  of  the  structure  are  linear  (linear  stress-strain  relationship).  The  slope  of  the 
curve,  3000  N/  kg,  was  calculated  from  the  tabulated  data  using  a  least  squares  fit. 

The  slope  of  the  curve  is  the  estimated  load  capability  per  kg  of  structural  mass.  As 
the  load  decreased  under  2(XX)0  N  the  minimum  gauge  restraint  influence  began  to 
dominate  the  minimum  weight  design.  The  minimum  gauge  constraint  was  enforced 
in  this  design  section.  With  the  minimum  gauge  restraint  active,  the  minimum  design 
weight  was  16.2  kg  (all  members  at  minimum  size).  As  long  as  the  design  load  is 
well  above  20000  N,  however,  a  good  guide  fm*  estimating  the  required  structural 
weight  for  a  given  design  load  is  to  divide  die  design  load  by  the  slope  (3000  N/kg)  of 
the  design  load/  design  weight  curve.  The  reader  is  reminded  that  this  curve  is  based 
on  rigid  loads  for  the  given  Mach  number/  angle  of  attack  pair.  The  effect  that 
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flexibility  would  have  on  the  optimization  process  could  not  be  determined  with  the 
capabilities  available.  It  would  be  interesting  to  know  how  the  flexible  design  curve 
would  compare  to  this  curve. 

The  next  relationship  that  was  explored  was  the  influence  that  design 
constraints  had  on  the  optimization  problem.  The  stress  constraint  case  provided  the 
lowest  objective  weight  (38.2  kg).  After  an  additional  minimum  gauge  constraint  was 
added  the  weight  increased  only  from  38.2  kg  to  38.4  kg.  One  of  the  reasons  for  such 
a  small  weight  penalty  is  the  ability  of  the  optimization  code  to  compensate  for  the 
size  limitation  by  reducing  the  sizes  of  other  elements  which  were  not  fully  stressed 
nor  at  the  minimum  gauge  value.  This  reduction  increased  the  forces  on  the  elements 
which  were  at  the  minimum  gauge,  therefore,  utilizing  the  potential  load  bearing 
capability  of  the  elements  at  a  minimum  gauge  value  (sec  figures  4-14  through  4-18). 

The  last  constraint  was  the  tip  deflection  constraint.  This  constraint  was  aimed 
at  reducing  the  expected  tip  acceleration.  This  could  not  be  achieved  by  merely 
constraining  the  displacement,  unfortunately,  because  the  square  of  the  first  fiequency 
of  vibration  increased  more  rapidly  than  the  tip  deflection  decreased.  (See  figure  S-2) 
This  is  not  a  complete  suiprise  since  fiequency  increases  with  additional  stiffness  (see 
equations  2-40  through  2-42)  and  displacement  decreases  with  additional  stifiness 
(equation  2-33).  If,  unlike  figure  5-2,  a  feasible  solution  does  exist  then  a  constraint 
formula  would  have  to  be  formulated  by  producing  a  sensitivity  equation  similar  to 
equations  2-33  and  2-42  except  the  sensitivity  equation  would  be  concerned  with 
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Figure  5*2:  Example  of  possible  displacement  and  natural  frequency  curves  as  a  function 
of  structural  stiffness. 


This  was  not  incoiporated  into  this  design.  The  difficulty  with 


this  type  of  constraint  is  the  requirement  to  utilize  both  the  strength  analysis  equations 
and  the  frequency  analysis  equations  (NASTRAN  can  not  perform  both  of  these 
simultaneously). 

The  displacement  constraint  increased  the  weight  of  the  design  by 
approximately  6  kg.  The  biggest  increases  in  design  variable  size  was  to  the  rib 
panels,  followed  by  smaller  increases  to  the  spar  and  skin  panels.  So  if  tip 
displacement  becomes  an  issue  in  the  tail  design,  the  structural  members  which  would 
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most  effect  the  displacement  are  the  rib  panels,  followed  by  the  skin  and  spar  panels. 


5.3  Aeroelastic  Analysis 


The  static  aeroelastic  response  of  the  three  constrained  designs  showed  that  the 
added  load  from  the  flexibility  affected  the  structural  response  less  than  three  percent 
except  with  the  stress  constrained  design  with  which  the  lift  coefficient,  , 

increased  over  seven  percent.  (See  Table  4-1)  In  each  case,  however,  the  lift 
coefficient  did  increase.  The  larger  increases  in  the  rolling  moment,  Cw  and  C»  , 

suggest  that  most  of  the  increase  in  lift  is  due  to  a  positive  rotation  near  the  tip  of  the 
tail.  There  were  no  stress  concentrations  in  any  of  the  elements  near  the  tip,  however, 
instead  there  was  a  general  increase  in  stress  values  throughout  the  members.  All 
member  stresses  remained  under  the  yield  stress  limit. 

A  dynamic  analysis  of  each  of  the  designs  showed  that  the  first  mode  was  a 
bending  mode  for  each  of  the  designs.  Only  in  the  stress  constrained  design  was  the 
second  mode  a  torsion  mode.  For  each  of  the  other  designs  the  second  mode  was  a 
bending  mode.  It  is  interesting  to  note  that  the  largest  aeroelastic  influence  came  from 
the  stress  constrained  design.  It  appears  that  the  modal  shapes  might  influence  the 
magnitude  to  which  the  aeroelastic  loads  influence  the  overall  loads. 

Two  relationships  were  determined  from  the  flutter  analysis;  the  variation  in 
flutter  speed  between  the  different  constraint  designs,  and  the  variation  in  flutter  speed 
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with  respect  to  mass  distribution.  The  largest  flutter  speed  for  the  three  design  cases 
was  with  the  minimum  gauge  case.  By  shifting  the  mass  distribution  of  the  minimum 
gauge  design  so  that  much  of  the  mass  was  in  the  forward  section  of  the  tail,  there 
was  a  large  gain  in  flutter  speed.  (500  m/s  — >  737  m/s)  By  shifting  the  mass  aft,  the 
flutter  speed  decreased  to  480  m/s.  The  flutter  speeds  are  based  on  the  aerodynamic 
influence  matrix  at  M=0.8.  None  of  the  flutter  speeds  are  in  the  defined  flight 
envelope  (M  <  0.8  sea  level  or  270  m/s).  An  important  fact  to  rotice  is  that  the  flutter 
phenomenon  appears  to  be  most  influenced  by  the  first  and  second  mode  of  each 
design.  These  modes  are  bending  modes,  except  in  the  stress  constraint  design  in 
which  the  second  mode  is  a  torsion  mode. 


5.4  Closing  Comments 

This  thesis  provides  a  base  for  additional  research  to  build.  The  analysis  of  the 
tail  utilized  structural  optimization  to  improve  the  design.  The  design  environment 
was  predicted  with  doublet-lattice  aerodynamics.  There  are  three  major  areas  fix)m 
which  the  tail  design  can  be  modified 

(1)  Structural  Model, 

(2)  Optimization  Parameters, 

(3)  Aerodynamic  theory. 

The  structural  model  can  be  modified  to  include  the  effect  of  topics  such  as  control 
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surfaces,  a  flexible  fuselage,  or  an  all-movable  tail  design.  How  would  the  aeroelastic 
response  vary  for  each  of  these  designs? 

The  area  with  the  most  attention  in  todays  literature  (35:982)  is  optimization. 
To  be  more  exact,  multi-disciplinary  optimization.  An  example  from  this  thesis  is  the 
inability  to  optimize  the  tail  members  with  respect  to  flexible  aerodynamic  loads.  The 
aerodynamic  force  equation  (Eqn  2-18)  needs  to  be  tied  into  the  influence  equations  of 
section  2-5.  Two  new  areas  which  are  gaining  a  lot  of  attention  in  the  engineering 
community  are  shape  sensitivities  where  the  influence  of  the  geometry  of  the  design  is 
incmporated  into  the  optimization,  and  the  ability  to  insert  or  remove  structural 
members  depending  on  whether  they  improve  or  degrade  the  objective  of  the 
optimization.  With  respect  to  aeroelasticity,  computer  code  has  been  developed 
which  allows  the  user  to  include  as  constraints,  or  an  objective,  stability  derivatives, 
flutter  speed,  divergence  speed,  etc.  More  and  more  relationships  are  being 
implemented  constantly.  As  researchers  and  programmers  incorporate  these 
relationships  into  their  computer  codes,  the  vertical  tail  model  could  be  used  as  a  test 
bed  for  evaluation. 

The  last  area  is  the  area  of  aerodynamic  modeling.  The  aerodynamic  model 

used  in  this  thesis  uses  linear  potential  flow  theory.  The  non-linear  influence  of 

high  angles  of  attack,  transonic  flow,  and  the  viscous  effects  of  turbulence  were 

ignored  in  the  current  aerodynamic  niodel.  The  influence  that  these  effects  might  be 

significant,  therefore,  a  new  approach  to  predicting  aerodynamic  loads  might  be 

warranted. 

\ 
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Appendix  A  :  Design  Modal  Shapes 


The  first  four  modal  shapes  of  the  five  designs  described  in  chapter  four  are 
given  in  the  following  pages.  Each  graph  displays  three  lines  which  depict  the  spars 
of  the  tail  box  design  as  seen  from  a  point  in  the  x-y  plane  looking  in  the  x-direction. 
The  front  most  spar  is  labeled  (FS).  the  mid-spar  is  labeled  (MS),  and  the  aft  spar  is 
labeled  (AS).  The  frequency  of  vibration  is  included  in  each  of  the  figure  titles  as 
"(Ci)=  ..)''.  These  drawings  should  illustrate  the  dominance  of  bending  or  torsion  in 
each  mode. 


A-1 


Figure  A<1:  First  mode  of  vibration  (to  ==  10.2  hz)  -  stress  constrained  case 
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Figure  A«2:  Second  mode  of  vibratitMi  (a)-41.3  hz)  -  stress  constrained  case 
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Nodal  Displacawnt  «  Displacement 
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AO:  Third  mode  of  vibration  (cd=S3.9  hz)  -  stress  constrained  case. 


Figure  A>4:  Fourth  mode  of  vibration  ((0=89  J  hz)  •  stress  constrained  case. 

AO 


Figure  A-5:  First  mode  of  vibration  ((0=11.3  hz)  -  minimum  gauge  case. 
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Rgure  A-6:  Second  mcxie  of  vibration  (to  =  43.01  hz)  -  minimum  gauge  case. 
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Hodal  Displacement  A  Modal  Displacement 
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A-7:  Third  mode  of  vibration  (a>=55.23  hz)  -  minimum  gauge  case. 
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Figure  A«8:  Fourth  mode  of  vibration  (0=97.7  hz)  •  minimum  gauge  case. 
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Modal  Displacement 


Figure  A-9:  First  mode  of  vibration  (tD=  16.91  hz)  -  displacement  constraint  case. 


Figure  A'lO:  Second  mode  of  vibration  (o)sS0.6  hz)  -  displacement  constraint  case. 
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Modal  DlaplaceMnt  <t  Displacement 


0.2 


\>11:  Third  mode  of  vibration  (a)=56.8  hz)  -  displacement  constraint  case. 
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Figure  A-12:  Fourth  mode  of  vibration  (ci)Bl06.3  bz)  •  displacement  constraint  case. 
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Figure  A-13;  Fi 
case. 
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Figure  A-14:  Se( 
case. 
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Figure  A-IS:  Thiid  mode  of  vibration  ((0=75.41  hz)  -  mass  forward  (minimum  gauge) 
case. 


Figure  A*16:  Fourth  mode  of  vibration  (oi^s!90.58  hz)  -  mass  forward  (minimum  gauge) 
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A-17:  First  mode  of  vibration  (a)=13.59  hz)  -  mass  aft  (minimum  gauge)  case. 
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Figure  A-18:  Second  mode  of  vibration  (o>s43.42  hz)  -  mass  aft  (minimum  gauge)  case. 


A-10 


Modal  Displacement  ^  Modal  Displacement 
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T'Coordinate  (n) 

Figure  A>20:  Fourth  mode  of  vibration  ((i)sl00.84  hz)  -  mass  aft  (minimum  gauge)  case. 
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Appendix  B  :  Input  Data  Decks 


Four  sample  data  decks  are  provided  to  demonstrate  the  procedure  and  to 
indicate  the  type  of  information  that  is  required  to  submit  the 

(a)  Structural  Optimization  Analysis, 

(b)  Static  Aeroelastic  Analysis, 

(c)  Flutter  Analysis. 

To  reduce  the  amount  of  repetition,  the  first  input  data  set  is  the  structural  defliution. 
This  definition  is  common  among  all  three  analysis.  The  lines  beginning  with  a  ’$’ 
are  comment  lines. 


Tail  Structural  Definition 


$  The  SPC  card  enforces  the  boundary  conditions  necessary  to  cantilever 
$  the  toot  of  the  tail. 

$ 

SPCl  1  123456  1  2  3  31  32  33 

$ 

$  The  MATl  card  defines  the  material  prcqierties  of  the  aluminum. 

$ 

MATl  17.24el0  0.33  2800.  +MAT 

-t-MAT  2.76E8  2.76E8  1.89E8 
$ 

$  Property  Cards  are  Identified  as  follows... 

$  PROD(l,23.4)  »  Posts 
$  PROD(S.6.7,8)  Spars 
$  PROD(9,10,11)  a  Ribs 
$  PSHELL(4)a  Skin 
$  PSHELL(S)«  SparShear 
$  PSHELL(Q>  RibShear 
$ 
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$  The  property  cards  give  the  cross-sectional  area  of  the  bar  elements 
$  and  the  thickness  of  panel  elements. 

$ 


PROD 

1 

1  0.01  1.0  1.0 

PROD 

2 

1  0.01  1.0  1.0 

PROD 

3 

1  0.01  1.0  1.0 

PROD 

4 

1  0.01  1.0  1.0 

PROD 

5 

1  0.01  1.0  1.0 

PROD 

6 

1  0.01  1.0  1.0 

PROD 

7 

1  0.01  1.0  1.0 

PROD 

8 

1  0.01  1.0  1.0 

PROD 

9 

1  0.01  1.0  1.0 

PROD 

10 

1  0.01  1.0  1.0 

PROD 

11 

1  0.01  1.0  1.0 

PSHELL 

12 

1  0.001  1  1. 

1 

0.0 

PSHELL 

13 

1  0.001  1  1. 

1 

0.0 

PSHELL 

14 

1  0.001  1  1. 

1 

0.0 

PSHELL 

15 

1  0.001  1  1. 

1 

0.0 

PSHELL 

16 

1  0.001  1  1. 

1 

0.0 

PSHELL 

17 

1  0.001  1  1. 

1 

0.0 

PSHELL 

18 

1  0.001  1  1. 

1 

0.0 

PSHELL 

19 

1  0.001  1  1, 

1 

0.0 

PSHELL 

20 

1  0.001  1  L 

1 

0.0 

PSHELL 

21 

1  0.001  1  1. 

1 

0.0 

GRDSET 

$ 

456 

$  GRID: 

Nodal  Point  Definition. 

$  CROD: 

Nodal  connectivity  of  the  bar  element. 

$  QUAD4: 

Nodal  Connectivity  of  the  Panel  element 

GRID 

1 

0  0.60000  0.000000.060000 

0 

GRID 

2 

0  1.10000  0.000000.080000 

0 

GRID 

3 

0  1.60000  0.000000.060000 

0 

GRID 

4 

0  1.29550  0244000.076747 

0 

GRID 

5 

0  1.01900  0.465400.055346 

0 

GRID 

6 

0  1 J0413  0.782230.052178 

0 

GRID 

7 

0  1.56129  0.576250.072309 

0 

GRID 

8 

0  1.84296  0.350610.056494 

0 

GRID 

9 

0  1.58925  1.099050.049009 

0 

GRID 

10 

0  1.82708  0.908500.067887 

0 

GRID 

11 

0  2.08593  0.701220.052988 

0 

GRID 

12 

0  1.87438  1.415880.045841 

0 

GRID 

13 

0  2.09286  1240750.063457 

0 

GRID 

14 

0  2.32889  1.051840.049482 

0 

GRID 

15 

0  2.15950  1.732700.042673 

0 

GRID 

16 

0  2.35865  U73000.0S9027 

0 

GRID 

17 

0  157185  1.402450.045976 

0 

GRID 

18 

0  2.44463  2.049530.039505 

0 
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GRID 

19 

0  2.62444  1.905250.054597 

0 

GRID 

20 

0  2.81481  1.753060.042469 

0 

GRID 

21 

0  2.72975  2.366350.036337 

0 

GRID 

22 

0  2.89023  2.237500.050167 

0 

GRID 

23 

0  3.05778  2.103670.038963 

0 

GRID 

24 

0  3.01488  2.683170.033168 

0 

GRID 

25 

0  3.15601  2.569750.045737 

0 

GRID 

26 

0  3.30074  2.454290.035471 

0 

GRID 

27 

0  3.30000  3.000000.030000 

0 

GRID 

28 

0  3.42180  2.902000.041307 

0 

GRID 

29 

0  3.54370  2.804900.031951 

0 

GRID 

31 

0  0.60000  0.00000-.060000 

0 

GRID 

32 

0  1.10000  0.00000-.080000 

0 

GRID 

33 

0  1.60000  0.00000-.060000 

0 

GRID 

34 

0  1.29550  0.24400-.076747 

0 

GRID 

35 

0  1.01900  0.46540-.055346 

0 

GRID 

36 

0  1.30413  0.78223-.052178 

0 

GRID 

37 

0  1.56129  0.57625-.072309 

0 

GRID 

38 

0  1.84296  0.35061..056494 

0 

GRID 

39 

0  1.58925  1.09905..049009 

0 

GRID 

40 

0  1.82708  0.90850..067887 

0 

GRID 

41 

0  2.08593  0.70122-.052988 

0 

GRID 

42 

0  1.87438  1.41588-.045841 

0 

GRID 

43 

0  2.09286  1.24075-.063457 

0 

GRID 

44 

0  2.32889  1.05184-.049482 

0 

GRID 

45 

0  2.15950  1.73270-.042673 

0 

GRID 

46 

02.35865  1.57300..059027 

0 

GRID 

47 

0  2.57185  1.40245-.045976 

0 

GRID 

48 

0  2.44463  2.04953-.039505 

0 

GRID 

49 

0  2,62444  1.90525-.054597 

0 

GRID 

50 

0  2.81481  1.75306-.042469 

0 

GRID 

51 

0  2.72975  236635-.036337 

0 

GRID 

52 

0  2.89023  2.23750..050167 

0 

GRID 

53 

0  3.05778  2.10367-.038963 

0 

GRID 

54 

0  3.01488  2.68317-.033168 

0 

GRID 

55 

0  3.15601  2J6975..045737 

0 

GRID 

56 

0  3.30074  2.45429-.035471 

0 

GRID 

57 

0  3.30000  3.00000-.030000 

0 

GRID 

58 

0  3.42180  2,90200-.041307 

0 

GRID 

59 

0  3.54370  2.80490..031951 

0 

CROD 

1 

5  1  5 

CROD 

2 

5  2  4 

CROD 

3 

5  5  6 

CROD 

4 

5  4  7 

CROD 

5 

5  3  8 

CROD 

6 

5  6  9 

CROD 

7 

5  7  10 

CROD 

8 

5  8  11 

CROD 

9 

6  9  12 

CROD 

10 

6  10  13 
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CROD 

11 

6 

CROD 

12 

6 

CROD 

13 

6 

CROD 

14 

6 

CROD 

15 

7 

CROD 

16 

7 

CROD 

17 

7 

CROD 

18 

7 

CROD 

19 

7 

CROD 

20 

7 

CROD 

21 

8 

CROD 

22 

8 

CROD 

23 

8 

CROD 

24 

8 

CROD 

25 

8 

CROD 

26 

8 

CROD 

27 

9 

CROD 

28 

9 

CROD 

29 

9 

CROD 

30 

9 

CROD 

31 

9 

CROD 

32 

9 

CROD 

33 

10 

CROD 

34 

10 

CROD 

35 

10 

CROD 

36 

10 

CROD 

37 

10 

CROD 

38  . 

10 

CROD 

39 

11 

CROD 

40 

11 

CROD 

41 

11 

CROD 

42 

11 

CROD 

43 

11 

CROD 

44 

11 

CTRIA3 

45 

12 

CQUAD4 

46 

12 

CQUAD4 

47 

12 

CQUAD4 

48 

12 

CQUAD4 

49 

12 

CQUAD4 

50 

12 

CQUAD4 

51 

13 

CQUAD4 

52 

13 

CQUAD4 

53 

13 

CQUA04 

54 

13 

CQUAD4 

55 

13 

CQUAD4 

56 

13 

CQUAD4 

57 

14 

CQUAD4 

58 

14 

CQUAD4 

59 

14 

CQUAD4 

60 

14 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 

4 

5 
7 

6 
10 
9 
13 
12 
16 
15 
19 
18 
22 
21 
25 
24 
28 
27 


4 

2 

0.0 

4 

5 

1 

0.0 

8 

7 

4 

0.0 

7 

6 

5 

0.0 

11 

10 

7 

0.0 

10 

9 

6 

0.0 

14 

13 

10 

0.0 

13 

12 

9 

0.0 

17 

16 

13 

0.0 

16 

15 

12 

0.0 

20 

19 

16 

0.0 

19 

18 

15 

0.0 

23 

22 

19 

0.0 

22 

21 

18 

0.0 

26 

25 

22 

0.0 

25 

24 

21 

0.0 
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11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

3 

4 
8 
7 
11 
10 
14 
13 
17 
16 
20 
19 
23 
22 
26 
25 
29 
28 

3 

2 

3 

4 
8 
7 
11 
10 
14 
13 
17 
16 
20 
19 
23 
22 


CQUAD4 

61 

14 

26 

29 

28 

25 

0.0 

CQUAD4 

62 

14 

25 

28 

27 

24 

0.0 

CROD 

101 

5 

31 

35 

CROD 

102 

5 

32 

34 

CROD 

103 

5 

35 

36 

CROD 

104 

5 

34 

37 

CROD 

105 

5 

33 

38 

CROD 

106 

5 

36 

39 

CROD 

107 

5 

37 

40 

CROD 

108 

5 

38 

41 

CROD 

109 

6 

39 

42 

CROD 

110 

6 

40 

43 

CROD 

111 

6 

41 

44 

CROD 

112 

6 

42 

45 

CROD 

113 

6 

43 

46 

CROD 

114 

6 

44 

47 

CROD 

115 

7 

45 

48 

CROD 

116 

7 

46 

49 

CROD 

117 

7 

47 

50 

CROD 

118 

7 

48 

51 

CROD 

119 

7 

49 

52 

CROD 

120 

7 

50 

53 

CROD 

121 

8 

51 

54 

CROD 

122 

8 

52 

55 

CROD 

123 

8 

53 

56 

CROD 

124 

8 

54 

57 

CROD 

125 

8 

55 

58 

CROD 

126 

8 

56 

59 

CROD 

127 

9 

33 

34 

CROD 

128 

9 

34 

35 

CROD 

129 

9 

38 

37 

CROD 

130 

9 

37 

36 

CROD 

131 

9 

41 

40 

CROD 

132 

9 

40 

39 

CROD 

133 

10 

44 

43 

CROD 

134 

10 

43 

42 

CROD 

135 

10 

47 

46 

CROD 

136 

10 

46 

45 

CROD 

137 

10 

50 

49 

CROD 

138 

10 

49 

48 

CROD 

139 

11 

53 

52 

CROD 

140 

11 

52 

51 

CROD 

141 

11 

56 

55 

CROD 

142 

11 

55 

54 

CROD 

143 

11 

59 

58 

CROD 

144 

11 

58 

57 

creiA3 

145 

12 

32 

34 

33 

0.0 

CQUAD4 

146 

12 

35 

34 

32 

31 

0.0 

CQUAD4 

147 

12 

37 

38 

33 

34 

0.0 

CQUAD4 

148 

12 

36 

37 

34 

35 

0.0 
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CQUAD4 

149 

12 

40 

41 

38 

37 

0.0 

CQUAD4 

150 

12 

39 

40 

37 

36 

0.0 

CQUAD4 

151 

13 

43 

44 

41 

40 

0.0 

CQUAD4 

152 

13 

42 

43 

40 

39 

0.0 

CQUAD4 

153 

13 

46 

47 

44 

43 

0.0 

CQUAD4 

154 

13 

45 

46 

43 

42 

0.0 

CQUAD4 

155 

13 

49 

50 

47 

46 

0.0 

CQUAD4 

156 

13 

48 

49 

46 

45 

0.0 

CQUAD4 

157 

14 

52 

53 

50 

49 

0.0 

CQUAD4 

158 

14 

51 

52 

49 

48 

0.0 

CQUAD4 

159 

14 

55 

56 

53 

52 

0.0 

CQUAD4 

160 

14 

54 

55 

52 

51 

0.0 

CQUAD4 

161 

14 

58 

59 

56 

55 

0.0 

CQUAD4 

162 

14 

57 

58 

55 

54 

0.0 

CROD 

201 

1 

35 

5 

caiOD 

202 

1 

34 

4 

CROD 

203 

1 

36 

6 

CROD 

204 

1 

37 

7 

CROD 

205 

1 

38 

8 

CROD 

206 

2 

39 

9 

CROD 

207 

2 

40 

10 

CROD 

208 

2 

41 

11 

CROD 

209 

2 

42 

12 

CROD 

210 

2 

43 

13 

CROD 

211 

2 

44 

14 

CROD 

212 

3 

45 

15 

CROD 

213 

3 

46 

16 

CROD 

214 

3 

47 

17 

CROD 

215 

3 

48 

18 

CROD 

216 

3 

49 

19 

CROD 

217 

3 

50 

20 

CROD 

218 

4 

51 

21 

CROD 

219 

4 

52 

22 

CROD 

220 

4 

53 

23 

CROD 

221 

4 

54 

24 

CROD 

222 

4 

55 

25 

CROD 

223 

4 

56 

26 

CROD 

224 

4 

57 

27 

CROD 

225 

4 

58 

28 

CROD 

226 

4 

59 

29 

CQUAD4 

227 

15 

31 

35 

5 

1 

0.0 

CQUAD4 

228 

15 

32 

34 

4 

2 

0.0 

CQUAD4 

229 

IS 

35 

36 

6 

5 

0.0 

CQUAD4 

230 

IS 

34 

37 

7 

4 

0.0 

CQUAD4 

231 

15 

33 

38 

8 

3 

0.0 

CQUAD4 

232 

15 

36 

39 

9 

6 

0.0 

CQUAD4 

233 

IS 

37 

40 

10 

7 

0.0 

CQUAD4 

234 

15 

38 

41 

11 

8 

0.0 

CQUAD4 

235 

16 

39 

42 

12 

9 

0.0 

CQUAD4 

236 

16 

40 

43 

13 

10 

0.0 
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CQUAD4 

237 

16 

41 

44 

14 

11 

0.0 

CQUAD4 

238 

16 

42 

45 

15 

12 

0.0 

CQUAD4 

239 

16 

43 

46 

16 

13 

0.0 

CQUAD4 

240 

16 

44 

47 

17 

14 

0.0 

CQUAD4 

241 

17 

45 

48 

18 

15 

0.0 

CQUAD4 

242 

17 

46 

49 

19 

16 

0.0 

CQUAD4 

243 

17 

47 

50 

20 

17 

0.0 

CQUAD4 

244 

17 

48 

5 

21 

18 

0.0 

CQUAD4 

245 

17 

49 

52 

22 

19 

0.0 

CQUAD4 

246 

17 

50 

53 

23 

20 

0.0 

CQUAD4 

247 

18 

51 

54 

24 

21 

0.0 

CQUAD4 

248 

18 

52 

55 

25 

22 

0.0 

CQUAD4 

249 

18 

53 

56 

26 

23 

0.0 

CQUAD4 

250 

18 

54 

57 

27 

24 

0.0 

CQUAD4 

251 

18 

55 

58 

28 

25 

0.0 

CQUAD4 

252 

18 

56 

59 

29 

26 

0.0 

CQUAD4 

253 

19 

35 

5 

4 

34 

0.0 

CQUAD4 

254 

19 

34 

4 

3 

33 

0.0 

CQUAD4 

255 

19 

36 

6 

7 

37 

0.0 

CQUAD4 

256 

19 

37 

7 

8 

38 

0.0 

CQUAD4 

257 

19 

39 

9 

10 

40 

0.0 

CQUAD4 

258 

19 

40 

10 

11 

41 

0.0 

CQUAD4 

259 

20 

42 

12 

13 

43 

0.0 

CQUAD4 

260 

20 

43 

13 

14 

44 

0.0 

CQUAD4 

261 

20 

45 

15 

16 

46 

0.0 

CQUAD4 

262 

20 

46 

16 

17 

47 

0.0 

CQUAD4 

263 

20 

48 

18 

19 

49 

0.0 

CQUAD4 

264 

20 

49 

19 

20 

50 

0.0 

CQUAD4 

265 

21 

51 

21 

22 

52 

0.0 

CQUAD4 

266 

21 

52 

22 

23 

53 

0.0 

CQUAD4 

267 

21 

54 

24 

25 

55 

0.0 

CQUAD4 

268 

21 

55 

25 

26 

56 

0.0 

CQUAD4 

269 

21 

57 

27 

28 

58 

0.0 

CQUAD4 

270 

21 

58 

28 

29 

59 

0.0 

ENDDATA 


Structural  Optimization  Analysis 


This  is  a  san^le  listing  from  the  21  design  variable  input  deck.  The  constraints  which 
are  active  are  stress,  minimum  gauge,  and  tip  deflection  (node  28).  The  applied  force 
is  the  rigid  aerodynamic  loads  for  Mach  number  =  0.4  and  the  angle  of  attack  =  30 
degrees. 


ID  OPTIMIZE,STATIC 

SOL  200  $  Static  Analyris  &  Optimization 
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TIME  120 
CEND 

TITLE  =  STRENGTH  OPTIMIZATION  OF  WINGBOX 
SPC=  1 

PARAM,  APPC,  STATICS 
LOAD  =31 
DISP  =ALL 
STRESS  =  ALL 
BEGIN  BULK 
$ 

$  MODEL  MEASURED  IN  SI  UNITS 

$  LENGTH  =  M 

$  FORCE  =  NEWTONS 

$  TIME  =  SECONDS 

$  GRAY  =  9.81  M/S^2 

$ 

$  DESr-N  MODEL 

$ 

$ 

$ 

$  DESVAR(U3.4)  =  Posts 
$  DESVAR(5.6.7.8)  =  Spars 

$  DESVAR(9.10.11)  =  Ribs 

$  DESVAR(12.13.14)  =  Skin 
$  DESVAR(15,16,17.18)  =  SparShear 
$  DESVAR(19;i0^1)  =  RibShear 

$ 

$  DESVAR  :  Design  Variable  Definition 
$  ID  Label  Start  Minimiun  Maximum 


$ 

$ 

Value  Value 

Value 

DESVAR 

1 

STll  .lOOE-04 

lE-5 

10. 

DESVAR 

2 

ST12  .151E-04 

l£-5 

10. 

DESVAR 

3 

ST13  .151E-04 

l£-5 

10. 

DESVAR 

4 

ST14  .851E-04 

l£-5 

10. 

DESVAR 

5 

ST15  .514E-04 

lE-5 

10. 

DESVAR 

6 

ST16  .731E-04 

l£-5 

10. 

DESVAR 

7 

ST17  .751E-04 

l£-5 

10. 

DESVAR 

8 

ST18  .598E-04 

l£-5 

10. 

DESVAR 

9 

ST19  .411E-04 

l£-5 

10. 

DESVAR 

10 

ST20  .801E-04 

l£-5 

10. 

DESVAR 

11 

ST21  .801E-04 

l£-5 

10. 

g 

DESVAR 

12 

SP22.227E-02 

lJE-3 

.1 

DESVAR 

13 

SP23  .lOOE-02 

l£-3 

.1 

DESVAR 

14 

SP24  .132E-02 

l£-3 

.1 

DESVAR 

IS 

SP2S  .161E-02 

l.E-3 

.1 

DESVAR 

16 

SP26  .19SE-02 

l£-3 

.1 

DESVAR 

17 

SP27  .169E-02 

l.E-3 

.1 
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DESVAR 

DESVAR 

DESVAR 

DESVAR 


18 

SP28  .lOlE-02 

l.E-3 

.1 

19 

SP29  .964E-02 

l.E-3 

.1 

20 

SP30  .992E-02 

l.E-3 

.1 

21 

SP31  .940E-02 

l.E-3 

.1 

$ 

$  DVPRELl  :  Corrolation  between  design  variables  and  structural 
$  properties  (Property  Cards) 

$ 


DVPRELl 

1 

PROD 

1 

4  lE-5 

+DP1 

+DP1 

1 

1.0 

DVPRELl 

2 

PROD 

2 

4  lE-5 

+DP2 

+DP2 

2 

1.0 

DVPRELl 

3 

PROD 

3 

4  lE-5 

+DP3 

+DP3 

3 

1.0 

DVPRELl 

4 

PROD 

4 

4  lE-5 

+DP4 

+DP4 

4 

1.0 

DVPRELl 

5 

PROD 

5 

4  lE-5 

+DP5 

+DP5 

5 

1.0 

DVPRELl 

6 

PROD 

6 

4  l£-5 

+DP6 

+DP6 

6 

1.0 

DVPRELl 

7 

PROD 

7 

4  lE-5 

+DP7 

+DP7 

7 

1.0 

DVPRELl 

8 

PROD 

8 

4  lE-5 

+DP8 

+DP8 

8 

1.0 

DVPRELl 

9 

PROD 

9 

4  lE-5 

+DP9 

+DP9 

9 

1.0 

DVPRELl 

10 

PROD 

10 

4  l£-5 

+DP10 

+DP10 

10 

1.0 

DVPRELl 

11 

PROD 

11 

4  l.E-5 

+DP11 

+DP11 

c 

11 

1.0 

DVPRELl 

12 

PSHELL 

12 

4  lE-3 

+DP12 

+DP12 

12 

1.0 

DVPRELl 

13 

PSHELL 

13 

4  l£-3 

+DP13 

+DP13 

13 

1.0 

DVPRELl 

14 

PSHELL 

14 

4  l£-3 

+DP14 

+DP14 

14 

1.0 

DVPRELl 

15 

PSHELL 

15 

4  l£-3 

+DP15 

+DP15 

15 

1.0 

DVPRELl 

16 

PSHELL 

16 

4  lE-3 

+DP16 

+DP16 

16 

1.0 

DVPRELl 

17 

PSHELL 

17 

4  lE-3 

+DP17 

+DP17 

17 

1.0 

DVPRELl 

18 

PSHELL 

18 

4  lE-3 

+DP18 

+DP18 

18 

1.0 

DVPRELl 

19 

PSHELL 

19 

4  lE-3 

+DP19 

+DP19 

19 

1.0 

DVPRELl 

20 

PSHELL 

20 

4  l£-3 

+DP20 

+DP20 

20 

1.0 

DVPRELl 

21 

PSHELL 

21 

4  1£.3 

+DP21 
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+DP21 

21 

1.0 

> 

DRESPl 

1 

SI 

STRESS 

PROD 

2 

1 

DRESPl 

2 

S2 

STRESS 

PROD 

2 

2 

DRESPl 

3 

S3 

STRESS 

PROD 

2 

3 

DRESPl 

4 

S4 

STRESS 

PROD 

2 

4 

DRESPl 

5 

S5 

STRESS 

PROD 

2 

5 

DRESPl 

6 

S6 

STRESS 

PROD 

2 

6 

DRESPl 

7 

S7 

STRESS 

PROD 

2 

7 

DRESPl 

8 

S8 

STRESS 

PROD 

2 

8 

DRESPl 

9 

S9 

STRESS 

PROD 

2 

9 

DRESPl 

10 

SIO 

STRESS 

PROD 

2 

10 

DRESPl 

11 

Sll 

STRESS 

PROD 

2 

11 

DRESPl 

12 

S12 

STRESS 

PSHELL 

7 

12 

DRESPl 

13 

S13 

STRESS 

PSHELL 

7 

13 

DRESPl 

14 

S14 

STRESS 

PSHELL 

7 

14 

DRESPl 

15 

S15 

STRESS 

PSHELL 

7 

15 

DRESPl 

16 

S16 

STRESS 

PSHELL 

7 

16 

DRESPl 

17 

S17 

STRESS 

PSHELL 

7 

17 

DRESPl 

18 

S18 

STRESS 

PSHELL 

7 

18 

DRESPl 

19 

S19 

STRESS 

PSHELL 

7 

19 

DRESPl 

20 

S20 

STRESS 

PSHELL 

7 

20 

DRESPl 

21 

S21 

STRESS 

PSHELL 

7 

21 

DRESPl 

22 

S22 

STRESS 

PSHELL 

8 

12 

DRESPl 

23 

S23 

STRESS 

PSHELL 

8 

13 

DRESPl 

24 

S24 

STRESS 

PSHELL 

8 

14 

DRESPl 

25 

S25 

STRESS 

PSHELL 

8 

15 

DRESPl 

26 

S26 

STRESS 

PSHELL 

8 

16 

DRESPl 

27 

S27 

STRESS 

PSHELL 

8 

17 

DRESPl 

28 

S28 

STRESS 

PSHELL 

8 

18 

DRESPl 

29 

S29 

STRESS 

PSHELL 

8 

19 

DRESPl 

30 

S30 

STRESS 

PSHELL 

8 

20 

DRESPl 

31 

S31 

STRESS 

PSHELL 

8 

21 

DRESPl 

32 

S32 

STRESS 

PSHELL 

9 

12 

DRESPl 

33 

S33 

STRESS 

PSHELL 

9 

13 

DRESPl 

34 

S34 

STRESS 

PSHELL 

9 

14 

DRESPl 

35 

S35 

STRESS 

PSHELL 

9 

15 

DRESPl 

36 

S36 

STRESS 

PSHELL 

9 

16 

DRESPl 

37 

S37 

STRESS 

PSHELL 

9 

17 

DRESPl 

38 

S38 

STRESS 

PSHELL 

9 

18 

DRESPl 

39 

S39 

STRESS 

PSHELL 

9 

19 

DRESPl 

40 

S40 

STRESS 

PSHELL 

9 

20 

DRESPl 

41 

S41 

STRESS 

PSHELL 

9 

21 

DRESPl 

43 

D43 

DISP 

3 

28 

DRESPl 

45 

W 

WEIGHT 

DESOBJ 

45 

w 

MIN 

$ 

$  Stress  Constraint  Definitiois 

$ 
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DCONSTR 

1 

ALL  -2.76E8 

2.76E8 

DCONSTR 

2 

ALL  -2.76E8 

2.76E8 

DCONSTR 

3 

ALL  -2.76E8 

2.76E8 

DCONSTR 

4 

ALL  -2.76E8 

2.76E8 

DCONSTR 

5 

ALL  -2.76E8 

2.76E8 

DCONSTR 

6 

ALL  -2,76E8 

2.76E8 

DCONSTR 

7 

ALL  -2.76E8 

2.76E8 

DCONSTR 

8 

ALL  -2,76E8 

2.76E8 

DCONSTR 

9 

ALL  -2.76E8 

2.76E8 

DCONSTR 

10 

ALL  -2.76E8 

2.76E8 

DCONSTR 

11 

ALL  -2.76E8 

2.76E8 

DCONSTR 

12 

ALL  -2.76E8 

2.76E8 

DCONSTR 

13 

ALL  -2.76E8 

2.76E8 

DCONSTR 

14 

ALL  -2.76E8 

2.76E8 

DCONSTR 

15 

ALL  -2.76E8 

2.76E8 

DCONSTR 

16 

ALL  -2.76E8 

2.76E8 

DCONSTR 

17 

ALL  -2.76E8 

2.76E8 

DCONSTR 

18 

ALL  -2.76E8 

2.76E8 

DCONSTR 

19 

ALL  -2.76E8 

2.76E8 

DCONSTR 

20 

ALL  -2.76E8 

2,76E8 

DCONSTR 

21 

ALL  -2.76E8 

2.76E8 

DCONSTR 

22 

ALL  -2.76E8 

2.76E8 

DCONSTR 

23 

ALL  -2,76E8 

2.76E8 

DCONSTR 

24 

ALL  -2.76E8 

2.76E8 

DCONSTR 

25 

ALL  -2,76E8 

2.76E8 

DCONSTR 

26 

ALL  -2,76E8 

2.76E8 

DCONSTR 

27 

ALL  -2.76E8 

2.76E8 

DCONSTR 

28 

ALL  -2.76E8 

2.76E8 

DCONSTR 

29 

ALL  .2.76E8 

2.76E8 

DCONSTR 

30 

ALL  -2.76E8 

2.76E8 

DCONSTR 

31 

ALL  -2.76E8 

2.76E8 

DCONSTR 

32 

ALL  -2.76E8 

2.76E8 

DCONSTR 

33 

ALL  -2.76E8 

2.76E8 

DCONSTR 

34 

ALL  -2.76E8 

2.76E8 

DCONSTR 

35 

ALL  -2.7ffi8 

2.76E8 

DCONSTR 

36 

ALL  -2.76E8 

2.76E8 

DCONSTR 

37 

ALL  -2.76E8 

2.76E8 

DCONSTR 

38 

ALL  -2.76E8 

2.76E8 

DCONSTR 

39 

ALL  -2.76E8 

2.76E8 

DCONSTR 

40 

ALL  -2.76E8 

2.76E8 

DCONSTR 

41 

ALL  -2.76E8 

2.76E8 

DCONSTR 

43 

ALL  '.10 

.10 

DSCREEN 

1  STRESS  -.1 

DSCREEN 

2 

DISP  -1. 

DOPIERM 

• 

2 

20 

$  •  • 

4> 

*  • 

* 

$  ANALYSIS  MODEL 
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$ 

$ 

PARAM  AUTOSPC  YES 
$ 

$  Nodal  Forces  from  aerodynamic  loads  at  M=0.4  aoa  =  10  deg.  This 
$  data  was  collected  from  aeroelastic  analysis. 

$ 


FORCE 

30 

5 

11502. 

0. 

0. 

1. 

FORCE 

30 

4 

13. 

0. 

0. 

1. 

FORCE 

30 

6 

2463. 

0. 

0. 

1. 

FORCE 

30 

7 

-3545. 

0. 

0. 

1. 

FORCE 

30 

8 

1752. 

0. 

0. 

1. 

FORCE 

30 

9 

5186. 

0. 

0. 

1. 

FORCE 

30 

10 

-1645. 

0. 

0. 

1. 

FORCE 

30 

11 

1159. 

0. 

0. 

1. 

FORCE 

30 

12 

4698. 

0. 

0. 

1. 

FORCE 

30 

13 

-1988. 

0. 

0. 

1. 

FORCE 

30 

14 

1223. 

0. 

0. 

1. 

FORCE 

30 

15 

4691. 

0. 

0. 

1. 

FORCE 

30 

16 

-1935. 

0. 

0. 

1. 

FORCE 

30 

17 

1177. 

0. 

0. 

1. 

FORCE 

30 

18 

4477. 

0. 

0. 

1. 

FORCE 

30 

19 

-1885. 

0. 

0. 

1. 

FORCE 

30 

20 

1071. 

0. 

0. 

1. 

FORCE 

30 

21 

4122. 

0. 

0. 

1. 

FORCE 

30 

22 

-1756. 

0. 

0. 

1. 

FORCE 

30 

23 

882. 

0. 

0. 

1. 

FORCE 

30 

24 

3467. 

0. 

0. 

1. 

FORCE 

30 

27 

1746. 

0. 

0. 

1. 

FORCE 

30 

25 

-1406. 

0. 

0. 

1. 

FORCE 

30 

26 

735. 

0. 

0. 

1. 

FORCE 

30 

28 

-810. 

0. 

0. 

1. 

FORCE 

30 

29 

-431. 

0. 

0. 

1. 

$  TOTAL  AERODYNAMIC  LIFT  =  34963.6  NEWTONS. 

*  *  *  *  *  *  0  0 

$  The  forces  ate  incremented  to  a  30  deg  aoa  by  multiplying  the  forces 
$  by  3.0  since  for  rigid  aerodynamics  the  fmce  is  linear  with  respect 
$  to  aoa. 

$  <...> 

LOAD  31  3.0  1.0  30 

$<***♦*♦*♦* 


Static  Aeroelasticity 

The  following  data  deck  is  a  sample  deck  from  the  aeroelastic  analysis  of  one 
of  the  designs  at  an  angle  of  attack  of  30  degrees  at  a  Mach  number  of  0.4.  (dynamic 
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pressure  of  11348  N/m^) 


NASTRAN 

ID  AERO.STATIC 

SOL  21  $  Static  Aero  Analysis 

TIMES 

COMPILE  SOL21  SOUIN=MSCSOU  NOLIST  NOREF 
$ 

$  This  alter  section  forces  the  output  to  split  the  aerodynamic 
$  loads  into  the  rigid  contribution  and  the  flexible  contribution. 

$  »»»» 

alter  1102  $  V66A  &  V67 

IF  (  N0RSET>-1  )  UPARTN  USET.KAAa/KaALL.,/AV’L7’R70  $ 

ALTER  1169  $  V66A 

MPYAD  KALX,  UX.  /ALOAD/  $ 

SMPYAD  MSLRJlT.TRX.ux,7ILOAD/4/////l  $ 

ADD  ILOAD,  ALOAD/  TLOAD/  $ 

MATGPR  GPL.  USET,  SIL,  ILOAD//’X7’L7  $ 

MATGPR  GPL,  USET.  SIL.  ALOAD/rXTLY  $ 

MATGPR  GPL,  USET.  SIL,  TLOAD//’X’/’L7  $ 

MATGPR  GPL.  USET.  SIL,  PLX//*XTL7  $ 
alter  1170$  V66A 
mpyad  kaall,  ul.  /aeload/  $ 

MATGPR  GPL,  USET.  Stt..  aeload//’X7’L7  $ 

$««««« 

CEND 

TITLE*  3  SPAR  SWEPT  WING 

SUBTI  *  STATIC  AERODYNAMICS  (Cantilevered  Wing) 

$ 

$ 

SEALL*  ALL 
SPC  *  1 

LOAD  =30  $  only  for  stab.  der.  (supported  structure  case) 

DISP  =  ALL 
STRESS  =  ALL 
AEROF  =ALL 
APRES  =ALL 
subcase  1 
TRIM  =  1001 
subcase  2 
trim  =  1(X)2 

$ 

BEGIN  BULK 

$  MOl^  MEASURED  IN  SI  UNITS 
$  LENGTH-M 

$  FORCE  >  NEWTONS 

$  TIME  -  SECONDS 

$  GRAY  -  9.81  M/S>'2 

$ 

$  STEADY  AERO  MODEL 
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$ 

$  for  sea  level  P  =  101325.  N/m'^2  dens  =  1.225  kg/m''3 
$  a  =  340.290  m/s 

$ 

PARAM  COUPMASS  1 
PARAM  GRDPNT  60 


Restraint/Lump  Mass  #############«######## 


SUPORT 

60 

345 

ASETl 

3 

4 

THRU  29 

ASETl 

3 

34 

THRU  59 

SPCl 

1 

126 

60 

SPCl 

1 

456 

4  thru 

29 

SPCl 

1 

456 

34  thru 

59 

RBE2 

301 

60 

123456  32 

31 

1  2 

3+RB2 

+RB2 

$ 

GRAY 

33 

10 

9.810  0.0 

0.0 

-1.0 

FORCE 

20 

60 

l.+ll 

0.0 

0.0  1.0 

LOAD 

30 

1.0 

1.0  10 

1.0 

20 

CONM2 

310 

60 

0  l.+ll 

+CM2 

+CM2 

l.+ll 

l.+ll 

PARAM  WTMASS.1019368 
PARAM  AUNTTS. 1019368 

Wing  Position  Axis  **••*♦***♦•*•*•**♦*****♦*♦**•*•*• 
$  CID  RID  A1  A2  A3  B1  B2  B3 
CORD2R  2  0.0  0.0  0.0  0.  0.  15.0  +CORD2 

$  Cl  C2  C3 
+CORD2  2.0  0.0  15.0 

$ 


$***•**••■ 

'**  Wing  Stability  Axis  *••♦♦••****•*****•**••••••**•*** 

$  CID 

RID 

A1  A2  A3  B1 

B2  B3 

$  * 

*  * 

*  •  ♦  • 

*  * 

CORD2R 

11  2 

1.1  0.0  0.0  1.1 

0.0  1.0  +CORD11 

$  Cl 

C2 

C3 

+CORD11 2.0  0.0 

1.0 

$ 

SYMXZ 

SYMXY 

AEROS 

11  2.4  6.0  5.40 

1 

CAEROl 

2001 

2000  2  10  8 

1+CARO 

+CARO 

0.0 

0.0  0.0  2.400  3.000  3.000  0.0  1.200 

PAEROl 

2000 

$ 

smmuummmmnmiiii  Trim  Conditions  mmmmniimttmmmHmi 

$ 


AESTAT 

500 

ANGLEA 

AESTAT 

SOI 

PITCH 

AESTAT 

503 

URDD3 

AESTAT 

504 

URDD5 

AESTAT 

SOS 

URDD4 

$ 
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TRIM  1001  0.40  11348.  ANGLEA  0.52360  PITCH  0.0 

$ 

$ 

$  5  deg  =  .08727  20  deg  =  .34910  35  deg  =  .61087 

$  10  deg  =  .17453  25  deg  =  .43633  40  deg  =  .69813 

$  15  deg  =  .26180  30  deg  =  .52360  45  deg  »  .78540 

$ 

$  AERODYNAMICS  =>  STRUCTURE 

$ 

SPUNEl  2000  2001  2001  2080  17  0. 

SETl  17  4  thru  29 

$ 

$ 

$  !!!!!!!  Insert  Structural  Data  Here. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 


Flutter  Analysis 

The  following  input  deck  provided  the  data  necessary  to  estimate  the  flutter 
speed  of  the  vertical  tail  design. 


ID  KMETHOD  JLUTTER 
SOL  75  $  FLUTIER  ANALYSIS 
TIME  10 
CEND 

TITLE  =  SWEPT  WING 

LABEL  =  FLUTTER  ANALYSIS  (K-Method) 

SEALL  =  ALL 
SPC=  1 
METHOD  =  2 

CMETHOD  =  4  $  K-Method 

FMETHOD  =  3 

SVEC  =  ALL  $  Vibration  Modes 

SDisp  s  all  $  Rutter  Modes 

$ 

BEGIN  BULK 

PARAM  COUPMASS  1 
PARAM  GRDPNT  60 
PARAM  LMODES  6 
$ 

j*******************  Wing  Position  Axis  **«'*************«'**4'*****<i"******<)' 
$  CID  RID  A1  A2  A3  B1  B2  B3 
CORD2R  2  0.0  0.0  1.0  0.  0.  15.0  +OORD2 

«  PI  po  pi 

+aMlD2  2.0  0.0  15.0 

$ 

$  UNSTEADY  AERO  CARDS 
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$ 


AERO 

0  1.  1.20  1.23  1  1 

CAEROl 

1001  1000  2  10  8 

1+CARO 

+CARO 

0.0  0.0  0.0  2.400  3.000  3.000 

0.0 

1.200 

PAEROl 

$ 

1000 

$  AERO-STRUCTURE 

$ 

SPUNEl 

8  1001  1001  1080  17  0. 

SETl 

17  4  Am  29 

$  FLUTTER  FLIGHT  CONDITION 

$ 

FLUTTER 

3  K  911  912  913  1 

4 

FLFACT 

911  1. 

$ 

Mach# 

FLFACT 

912  .8 

$ 

$  for  varying  MACH  #  change  (1)  MKAEROl  and  (2)  FLFACT  912 


$ 

FLFACT  913  0.99 

.800  .700 

.600 

.500 

.450 

.400+FLFl 

+FLF1  .360  .340 

.320  .300 

290  .280 

.270 

.260+FLF2 

+FLF2  .255  .250 

.245  .240 

235  230 

.225 

220+FLF3 

+FLF3  .215  .205 

.190  .180 

.170  .160 

.150 

.140 

MKAEROl  .8 

+MKAROa 

+MKAROa  1.00  .700 

.600  .500 

.400 

.320 

.300  .290 

MKAEROl  .8 

+MKAROb 

+MKAROb  280  270 

.260  .240 

.225 

.210 

.205  200 

MKAEROl  .8 

+MKAROC 

+MKAROC  .195  .190 

.180  .170 

.160 

.150 

.140  .130 

EIGR  2  GIV 

EIGC  4  HESS 

+EGC 

MAX 

6 

6 

+EGC 

SPCl  1  123 

1  2 

3  31 

32 

33 

SPCl  1  123 

60 

$  restrict  modal  analysis  to  displacements  in  the  *3’  coordinate  direction. 
$ 

ASETl  3  4  THRU  29 

ASETl  3  34  THRU  59 

$ 

$  Concentrated  masses  at  selected  nodal  points. 

$ 


CONM2 

501 

5 

2. 

0.0 

CONM2 

502 

6 

2. 

0.0 

CONM2 

503 

9 

2. 

0.0 

CONM2 

504 

12 

2. 

0.0 

CQNM2 

SOS 

15 

2. 

0.0 

CONMl 

506 

18 

2. 

0.0 

CONM2 

507 

21 

2. 

0.0 

B-16 


CONM2 

508 

24 

2. 

0.0 

CONM2 

509 

27 

2, 

0.0 

CONM2 

510 

36 

2. 

0.0 

CONM2 

511 

42 

2. 

0.0 

CONM2 

512 

48 

2. 

0.0 

CONM2 

513 

54 

2. 

0.0 

CONM2 

514 

10 

2. 

0.0 

CONM2 

515 

16 

2. 

0.0 

CONM2 

516 

22 

2. 

0.0 

CONM2 

517 

29 

2. 

0.0 

$ 

$ 

$!!!!!!!!!!!!  Put  Structural  Definition  Data  in  here.!!!!!!!!!!!!!*!!! 
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